Accepted to ApJ: October 29, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 03/07/07 



ATMOSPHERIC ESCAPE FROM HOT JUPITERS 

Ruth A. Murray-Clay^'^, Eugene I. Chiang\ & Norman Murray^'"* 
Accepted to ApJ: October 29, 2008 

ABSTRACT 

Photoionization heating from UV radiation incident on the atmospheres of hot Jupiters may drive 
planetary mass loss. Observations of stellar Lyman-a absorption have suggested that the hot Jupiter 
HD 209458b is losing atomic hydrogen. We construct a model of escape that includes realistic heating 
and cooling, ionization balance, tidal gravity, and pressure confinement by the host star wind. We 
show that mass loss takes the form of a hydrodynamic ("Parker") wind, emitted from the planet's 
dayside during lulls in the stellar wind. When dayside winds are suppressed by the confining action of 
the stellar wind, nightside winds might pick up if there is sufficient horizontal transport of heat. A hot 
Jupiter loses mass at maximum rates of ~2 x 10^^ gs~^ during its host star's pre-main-sequence phase 
and ~2 x 10^°gs~^ during the star's main sequence lifetime, for total maximum losses of ~0.06% and 
~0.6% of the planet's mass, respectively. For UV fluxes i^uv ^ lO'^ergcm"^ s~^, the mass loss rate 
is approximately energy- limited and scales as M oc -Fxjv- For larger UV fluxes, such as those typical 
of T Tauri stars, radiative losses and plasma recombination force M to increase more slowly as F^^- 
Dayside winds are quenched during the T Tauri phase because of confinement by overwhelming stellar 
wind pressure. During this early stage, nightside winds can still blow if the planet resides outside the 
stellar Alfven radius; otherwise, even nightside winds are stifled by stellar magnetic pressure, and mass 
loss is restricted to polar regions. We conclude that while UV radiation can indeed drive winds from hot 
Jupiters, such winds cannot significantly alter planetary masses during any evolutionary stage. They 
can, however, produce observable signatures. Candidates for explaining why the Lyman-a photons of 
HD 209458 are absorbed at Doppler-shifted velocities of ±100 km/s include charge-exchange in the 
shock between the planetary and stellar winds. 

Subject headings: planetary systems — hydrodynamics — stars: individual (HD 209458) 



1. INTRODUCTION 

About 1/5 of the approximately 200 extrasolar planets 
discovered to date have masses comparable to Jupiter's, 
but orbit their host stars at distances less than 0.1 AU 
(e.g., [Butler et al, 2006). Since stellar heating inhibits 
the formation of g as giants aro und Sun-like stars at 
such distances fe.g.. lRafikovll2006l ). these "hot Jupiters" 
likely migrated i nward by disk torques f rom where they 
were born (e.g., iPapaloizou et al.l l2007l and references 
therein). 

Once parked (possibly because the disk wa s trun - 
cated by the stellar magnetosphere; iLin et al] I1996D . 
hot Jupiters are bathed in ultraviolet (UV) radiation 
from their host stars. Atmospheric gas is heated 
by photoionization of hydrogen, and escapes. Sev- 
eral groups have argued that the outfiows can evapo- 
rate g as giants nearly entirely, laying bare thei r rock y 
cores ()Lammer et al.l l2003l : iBaraffe et all l2004 I2005D . 
In fact, hot Jupiters are observed to have systemati- 
cally lower masses than extrasolar planets at la rger dis- 
tances fro m their st a rs ([ Zu cker fc Mazehl l2002( ) . How- 
ever, Hubbard et al. I (|2007aUbl l are unable to reproduce 
the mass distribution of hot Jupiters by experimentation 
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with mass loss histories. 

Observations of HD 209458b, th e first hot Jupiter 
observed to transit its host star (jHenrv et al.l l2000l : 
ICharbonneau et al]l2000l ). suggest t hat the planet may 
indeed be losing atomic hydrogen. IVidal-Madiar et aLl 
l|2003f ) used the Hubble Space Telescope Imaging Spec- 
trograph (STIS) in a high spectral resolution mode to 
measure Lyman-a emission from HD 209458b's host star 
in and out of transit. They observed, with 3cr confi- 
dence, the Lya flux — at wavelengths shifted from line 
center by Doppler equivalent velocities of ±100kms~^ — 
to decrease during transit by ~15% (as integrated over 
Doppler equivalent velocities extending approximately 
from -50 to -140 km/s and from -1-30 to -1-100 km/s). The 
authors attributed this decrease to absorption by inter- 
vening atomic hydrogen surrounding the planet: a halo 
of gas sufficiently distended and traveling at fast enough 
velocities to no longer be bound to the planet. The ab- 
sorption signal was reported to be stronger at blueshifted 
wavelengths than at redshifted wavelengths; blueshifted 
velocities (towards the observer, away from the star) were 
argued to arise from stellar radiation pressure acting on 
gas via the L yct lin e. 

Ben- Jaffel pOOTl ) reanalyzed the STIS data and agreed 
that stellar Lya photons were absorbed by planetary gas 
during transit, but found the signal to be weaker: the 
fiux decreased by only 8.9 ± 2.1%, at similar Doppler 
equivalent velocities of ±100 km/s. Moreove r, no pref- 
erence for blueshifted absorption was found. iBen- Jaffel 
(|2007( ) cautioned that intrinsic stellar variability could 
easily produce a spurious preference for either negative 
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or po sitive veloc i ties. Despite this and other comphca- 
tions. lBen-Jaffell (|2007( ) decided nonetheless that the ab- 
sorption signal could be of planetary origin. From the 
effective occulting area corresponding to an 8.9% drop 
in flux, it was concluded that the obscuring hydrogen, 
while occupying a "corona" significantly more inflated 
than the visible photosphere, remained bound within the 
pl anet's Roche lobe. This las t argument, as pointed out 
bv lVidal-Madiar et all (|2008D . is not secure. If the wave- 
length shifts associated with the putative absorption re- 
flect Doppler shifts from bulk flows, then the relevant 
velocities, regardless of whether they are positive or neg- 
ative, are larger than the planet's escape velocity and 
cannot arise from a bound hydrogen atmosphere. Fur- 
thermore, as a point of principle, gas can flow past the 
Roche lobe of the planet and elude detection if it is op- 
tically thin to Lya photons; therefore arguments based 
on effective occulting area are not con clusive. 

Certainly we share the concern of IBen-Jaffell ([20071 ) 
that stellar variability can corrupt any interpretation of a 
planetary wind. While .Vidal-Madiar et a l. (2003, 200^) 
attempted to account for this statistically, the out-of- 
transit time baseline may be too short to characterize 
confidently the intrinsic variability of the stellar Lya line 
(J. Winn, personal communication). There are further 
problems. Lyman-a absorption is measurable only in 
the line wings because of confusion near line center from 
the interstellar medium and from geocoronal (terrestrial) 
emission. Even the line wings, however, can be contam- 
inated by high ly time - variable geocoronal emission. 

Vidal-Madiar.eLalJ (f2004) further used STIS in a 
lower spectral resolution mode to measure wavelength- 
integrated fluxes in various lines. These measurements 
enjoyed greater sensitivity to variations in a spectral line 
at the cost of not resolving the line profile. These au- 
thors observed, with 2-3cr confidence, the H I Lya line 
flux to decrease by ~5% in transit. Purely from an ob- 
servational standpoint, this spectrally unresolved mea- 
surement is claimed to be consistent with the sp e ctrall y 
resolved measurement (IVidal-Madjar et all [20031 l2004f ). 
Despite the claimed consistency, the two sets of mea- 
surements are made more than two years apart, which 
is perhaps worrisome given how stellar chromospheric 
emission is highly time variable. Moreover, it is, of 
course, impossible to verify whether the spectrally unre- 
solved observations pertain to the same strikingly large 
Doppler equivalent velocities of -I-/- 100 km/s that are so 
clearly implicated in the spectrally resolved observations. 
In any case, decrements in O I and C II line emission 
were also observed, with similarly marginal confidence. 
Unfortunately, none of these tantalizing measurements 
could be reproduced after the STIS instrument failed in 
2004. lEhrenreich et al.l ()2008D recently attempted simi- 
larly spectrally unresolved measurements with the Hub- 
ble Advanced Camera for Surveys (ACS). Their detection 
of H I absorption in transit is uncertain but consistent 
with previous claims. 

Even if the STIS observations do signify a plane- 
tary outflow, the mass loss rate implied is not neces- 
sarily lar ge enough to seriou s ly red uce the mass of the 
planet. I Vidal-Madiar et all ()2003f ) claim a mass loss 
rate of 10^°gs~^ based on a fit to the absorption depth 
at large Dopper equivalent velocities vs. time, and on 
considerations of how radiation pressure can shape the 



cloud of hydrogen as it expands away from the planet. 
Note that radiation pressure is not necessarily claimed 
by these authors to drive the outflow; they acknowl- 
edge the need for hydrogen to be heated by the star 
and to expand to di stances approaching the Roche lobe 
(jVidal-Madjar fc Etang s 2004).^ Theoretical models of 
thermal winds heated by photoionization generally pro- 
duce mass loss rates that are several times 10^°gs~ ^ 
(lYellell2004[20MlTian et~an[2005l : lCarcia Mufiozl[200l . 
Over the several Gyr age of the system, somewhat less 
than 1% of a hot Jupite r's mass would be car ried away 
by such thermal winds. iGarcia Munod (|2007t l presents 
particularly convincing hydrodynamic escape models. 

Nevertheless, the question remains whether the vari- 
ous STIS observations of HD 209458b are correctly ex- 
plained as a planetary wind driven by photoionization. 
It is also unclear whether hot Jupiters lose significant 
mass early in their evolution, when the UV luminosi- 
ties of their host stars are enormously higher than their 
main-sequence values. 

In this paper, we demonstrate from a first-principles 
calculation that a hot Jupiter cannot lose a significant 
fraction of its mass via a planetary outfiow driven by UV 
photoionization at any stage during its lifetime, includ- 
ing the pre-main-sequence phase. We further show that 
the spectrall y resolved Lyman-a trans i t observations o f 
HD 209458b (| Vidal-Madiar et al.l[2003l : iBen Jaffe]|[200l 
probe velocities too large to refiect the bulk flow of a 
hot Jupiter wind. Observations of Lyman-a absorption 
at Doppler-equivalent velocities of ±100 km/s require ei- 
ther an additional sour ce of high- velocity H atoms (e.g., 
iHolmstrom et al.ll2008D or an enhancement in the den- 
sity of neutrals in or around the planetary wind (e.g., 
lBen-Jaffell2008[ ). We discuss these possibflities in ^ 

Our standard model is essentially that of a thermal or 
"Parker" wind — a flow accelerated by gas pressure from 
subso nic to superso nic velocities through a critical sonic 
point l| Parked [195 81 ) — with the added complication that 
the heating is external, from stellar UV irradiation. Our 
model includes realistic heating and cooling, ionization 
balance, and tidal gravity, but is simple enough that we 
can elucidate the basic physics and write down approxi- 
mate analytic formulae for the mass loss rate. 

[Garcia Munod (|2007l ) notes correctly that the plane- 
tary outflow need not take the form of a transonic wind.^ 
He points out that the host star wind can pressure- 
confine the planetary outfiow down to a subsonic breeze. 
Stellar wind interactions furt her complicate the p lane- 
tary fiow at large distances (cf. ISchneiter et al.|[2007l who 
model stellar wind interactions with an entirely neutral 
wind escaping from HD 209458b at large velocities). By 
using our hydrodynamics code to compute breeze solu- 
tions, we characterize qualitatively the extent to which 
outfiows from hot Jupiters are suppressed, both in the 
case of a wind emitted by a main-sequence Sun-like star, 

^ In ^3.51 we show that radiation pressure acting on hydrogen 
through the Lya transition adds at most 1% to the maximum mass 
loss rates achievable from a thermal wind heated by photoioniza- 
tion. 

^ In our terminology, "transonic" refers to the Parker wind which 
transitions fro m subsonic to mildly supersonic velocities, while 
IGarcia Munozi 1)200 7') uses the word to describe a breeze whose 
peak velocity is nearly sonic but which eventually decelerates to 
zero velocity at infinite distance. For the wind he reserves the 
word "supersonic." 
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and in the case of a pre-main-sequence T Tauri wind. 

In !j2l after reviewing some of the basic orders of magni- 
tude characterizing hot Jupiter atmospheres, we present 
our standard model of a steady transonic wind. We solve 
the equations of ionization balance and of mass, momen- 
tum, and energy conservation using a relaxation code, 
and we demonstrate the robustness of our solution by 
exploring a variety of input boundary conditions. The re- 
sults of our model are presented for UV fluxes spanning 
four orders of magnitude, ranging from those expected 
for quiet main-sequence solar analogs to those emitted by 
active T Tauri stars. Section [3] contains a wide-ranging 
discussion of various aspects of our wind solution. Of es- 
pecial interest are how the physics of mass loss changes 
as the UV flux increases from low main-sequence values 
to high pre-main-sequence values f H3.2p : how host star 
winds can squash dayside planetary outflows, thereby 
perhaps energizing nightside outflows (i )3.4p : and how 
radiation pressure is ineffective compared to UV pho- 
toionization heating in driving an outflow (i i3.5p . Finally, 
2] summa rizes our find ings, pinpoints why the estimate 
of Lamme r et al.l (|2003l) an d subsequent deter minations 
bv .Baraffe et al.l (|2004f ) and lBaraffe erall (p0051 of mass 
loss rates are erroneously high, and assesses the STIS 
observations, both spectrally resolved and unresolved. 

2. THE MODEL 

We construct a simple ID model for photoevapora- 
tive mass loss from a hot Jupiter. We focus on the flow 
originating from the substellar point on the planet, and 
assume that mass loss occurs in the form of a steady, 
hydrodynamic, transonic wind. These restrictions im- 
ply that we will calculate a maximum flux of mass from 
the planet, insofar as (a) the substellar point receives 
the maximum UV flux from the star, (b) tidal grav- 
ity acts most strongly along the substellar ray to accel- 
erate gas away from the planet, and (c) the transonic 
wind carries more mass than pressure-conflned, subsonic 
breezes.^ Applying our solution for the substellar, tran- 
sonic streamline over all Att steradians yields a hard up- 
per limit on the total rate of photoevaporative mass loss. 
How closely the actual rate of mass loss approaches this 
upper limit is discussed in Sj3l 

We assume the base of the wind is composed of atomic 
hydrogen and calculate how the flow becomes increas- 
ingly ionized as it approaches the star. We neglect the 
molecular chemistry of hydrogen and do not capture the 
H2/H dissociation front. This simpliflcation can be jus- 
tifled by showing that the temperature of the wind is 
higher than the ^2000 K required to thermally dissoci- 
ate H2, and by showing that above the r = 1 surface to 
photoionization, our solution is insensitive to our chosen 
boundary conditions. This wc do in Appendix [Al a sum- 
mary is given at the end of H2.'2.'2\ We further neglect 
helium and metals. We comment on the implications of 
this omission in S|4l 

The rest of this section is organized as follows. In ^2.1\ 
we write down the basic steady-state equations of mass, 
momentum, energy, and ionization balance. The numer- 
ical methods used to solve these coupled ordinary dif- 
ferential equations are detailed in ^2.21 which includes a 

For an introduction to Parker's II1958I) theory of transonic 
winds and subsonic bre ezes, see, e.g., "Introducti on to Stellar 
Winds," the textbook bv lLamers &: Cassinellil I I1999I ). 



listing of our boundary conditions. In i )2.3l we present 
the results of the model. For those wishing to skip to the 
punchline, a simple analytical description of our results 
may be found in 

To help orient the reader, we now supply some of our 
standard model parameters, together with several order- 
of-magnitude estimates characterizing the wind. At Ly- 
man continuum wavelengths, the solar UV luminosity is 
roughl y 10~^Lf7i, where L p^ is the bolometric solar lumi- 
nosity ijWoods et al.lll998[ ). To the extent that host stars 
of hot Jupiters are like the Sun, a hot Jupiter with an or- 
bital semi-major axis of a = 0.05 AU receives a UV flux of 
^uv — 450ergcm~^ between photon energies of 13.6 
eV and 4 eV (this is nearly identical to the flux em- 
ployed bv lGarcia Munoj ()2007D but we derive ours inde- 
pendently of that study). T his flux characterize s "mod- 
erate to low solar activity" (jWoods et al.lll998l ). and is 
not averaged over the planetary surface; for a discussion 
of the effects of surface averaging, see §3.31 We take the 
planet to have mass Mp = 0.7 Mj = 10'^° g and a flducial 
1-bar radius Rp = 1.4i?j — 10^° cm, where Mj and Rj 
are, respectively, the mass and radius of Jupiter. The 
planet's surface gravity g 700 cm/s^ is approximately 
the same as g ~ 10^ cm/s^ on Earth. We take the effec- 
tive radius of the planet's Roche lobe, inside of which the 
planet's gravity dominates the host star's tidal gravity, to 
equal the approximate distance to the planet's LI point: 
i?Rocho = [Mp/(3M*)]1/3q ^ 4.5i?p, where = Mq is 
the mass of the star. In many astrophysical situations — 
including ours, as will be shown — photoionized gas cools 
by radiation from coUisionally excited atomic hydrogen, 
which thermostats the gas temperature T to '^10^ K. 
The corresponding sound speed is '--^lO km/s. The hydro- 
static pressure scale height is H = kT/{m^g) ^ 0.1i?p, 
where k is the Boltzmann constant and mn is the mass 
of the hydrogen atom (of course the wind is not strictly 
hydrostatic but it is nearly so near its base where speeds 
are still subsonic). To travel a distance Rp at the sound 
speed takes a few hours. 

Finally, we can estimate the gas density and pressure 
where the wind is launched, i.e., where the bulk of the 
stellar UV radiation is absorbed. At a photon energy 
of hiyQ = 20 eV, the cross section for photoionization of 
hydrogen is ai,^ — 6 x 10~^^(/iJ/o/13.6 eV)~^ cm^ (e.g., 
ISpitzed [l978l ): optical depth unity is achieved in a neu- 
tral column A^H = 1/fi/o = 5 x lO^^cm"^; dividing this 
column by the scale height H gives a neutral density 
no 6 X 10® cm~^ (equivalently, a neutral mass density 
p ^ 10^^^ gcm^'^); and multiplying by kT gives a partial 
pressure P ~ 1 nanobar at the base of the wind. By con- 
trast, visible radiation from the star is absorbed at pres- 
sures closer to 1 bar, setting the temperature below the 
base of the wind to be Tbeiow ~ lO'^ K and the pressure 
scale height to be Hbeiow = fc2beiow/(2m.H5') 0.005i?p, 
where the factor of 2 accounts for the fact that the hy- 
drogen at depth is molecular. The smallness of i?boiow 
means that the wind is launched at a radius very nearly 
equal to Rp : reducing the pressure from 1 bar to 1 nano- 
bar takes about 20 scale heights or 0.1i?p. In other words, 
the radius at which UV photons are absorbed is approx- 
imately l.li?p. This radius enters signiflcantly into the 
magnitude of the mass loss rate, as we discuss in H3.2\ 
fgl and Appendix El 
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2.1. Basic Equations 

As stated above, we concentrate on the streamline orig- 
inating from the substellar point on the planet. From 
mass continuity, 



d_ 

dr 



(r^pv) = , 



(1) 



where r is the distance from the center of the planet to 
the star, and the gas has density p and velocity v. In 
the frame rotating with the planet's orbital frequency, 
momentum conservation implies 



dv 
or 



dP 
dr 



GMpp iGA'hpr 



(2) 



where G is the gravitational constant. We call the last 
term on the right-hand side of ([2]) the "tidal gravity" 
term: it is the sum of the centrifugal force and differential 
stellar gravity, along the ray joining the planet to the star 
(neglecting the sm all shift in the system barycenter away 
from the star; cf. iGarcfa Munoj|2007[ ). For simplicity, 
we neglect the Coriolis force, the magnitude of which is 
comparable to the magnitude of other forces only at the 
outer boundary of our calculation, near the Roche lobe 
radius where gas moves at approximately the planet's 
escape velocity. 
The equation for energy conservation is 



d 
or 



kT 



(7 - 1)/^ 



kTv dp 
p dr 



r + A 



(3) 



where the left-hand side tracks changes in the internal 
thermal energy of the fluid, p. is the mean molecular 
weight, and 7 = 5/3 is the usual ratio of specific heats 
for a monatomic ideal gas. On the right-hand side we 
have three terms denoting, respectively, cooling due to 
PdV work done by expanding gas, heating from pho- 
toionization, and cooling from radiation and conduction. 
We do not include a term proportional to the chemi- 
cal potential because changes in energy due to changes 
in the number of particles are already accounted for in 
our photoionization term F. Equation Q follows in a 
straightforward way from the standard steady-state en- 
ergy equation, V • (puv) — — PV • v -I- F -I- A, after using 
pvr^ = constant and u = kT/[{j — l)p] for the specific 
internal energy. 

We assume for simplicity that the UV flux is concen- 
trated at one photon energy hi^o — 20 eV. Then 



(4) 



where no is the number density of neutral H atoms, 
e — [huQ — 13.6eV)//iJ^o is the fraction of photon energy 
deposited as heat,* and 



uq dr 



(5) 



is the optical depth to ionizing photons. 

Equation (U) assumes that photoelectrons share their 
kinetic energy with other gas species locally. We have 
verified a posteriori that this assumption is justified. For 
our standard model, at the wind base, a photoelectron 

® This is a maximum efRciency and can be somewhat smaller if 
there are other ways for the primary ph otoelectron to dep osit its 
energy (e.g., by secondary ionization; cf . IWaite et al.|[l983l) . 



travels A,nfp ~ 2 x 10 ^R-p before colliding with a neutral 
H atom, assuming a cross section of ~10~^^cm^. The 
number of collisions required for the electron to give up 
most of its energy to surrounding H atoms is '^mn/nic, 
where me is the electron mass; the corresponding dis- 
tance random walked is (mH/mo)^/^Ainfp 7 x 10~'^i?p. 
A similarly short distance obtains at the outer periphery 
of our calculation, near the Roche lobe, where photo- 
electrons share their energy with H''" ions via Coulomb 
collisions. The margin of safety is still larger for our high 
flux model ( §2.3.21) . for which ion densities are greater. 

The main contribution to A is Lya radiation, emit- 
ted by neutral H atoms that are coUisionally excited by 
electrons: 

A « Alv. = -7.5 X 10-i9n+noe^"''''''/^ergcm-3s-i 

(6) 

where n+ is the number density of {— the num- 
ber densit y of electroii s), and all densities are measured 
in cm~^ (jBlackl [l981|) . Other cooling mechanisms — 
collisional ionization, radiative recombination, free-free 
emission, and thermal conduction — are negligible, as 
shown in ij2.3l (see also Appendix [B| . Our assumption 
that Lya photons are able to escape and thereby cool 
the flow is validated in Appendix ICl 

In ionization equilibrium, the rate of photoionizations 
balances the rate of radiative recombinations plus the 
rate at which ions are advected away: 



1 d 



no — r (T„g = n^arcc + -^-^ {r^n+v) , (7) 

hvo r'^ dr 

where arec = 2.7 x 10-"(r/10'' K)-0-9 is the Case 
B radiative recombinatio n coefficient for hydrogen ions 
(jStorev fc Hummei1ll995f ). By continuity, the advection 
term can be rewritten as 



(8) 



r dr dr 

where the ionization fraction /+ = nj^jn and n = ?i++rto 
is the total number density of hydrogen nuclei. Note that 
p — toh/(1 + /+). Collisional ionization is negligible 
compared to photoionization, as demonstrated in ii2.3l 

2.2. Numerical Method 

The problem of finding the structure of the wind is a 
two point boundary value problem (for an introd uction 
into the nature of such problems, see, e.g., iPress et al.l 
Il992f ). The two points are the base of the flow and the 
sonic point. We have solved the problem by construct- 
ing a relaxation code. In our case relaxation methods 
are preferred over shooting methods because for every 
transonic wind s olution there are an infinite number of 
breeze solutions (|Parkeij[l958f) . Furthermore, the sonic 
point is a critical point where derivatives, if not care- 
fully computed, can become singular. Instead of search- 
ing exhaustively in a multidimensional space for the one 
solution that "threads the needle" of the critical sonic 
point, it is more efficient to start with an approximate 
solution that already satisfies the sonic point conditions, 
and refine that solution to higher accuracy. Previous 
attempts that did not use relaxation algorithms to find 
transonic winds found, n ot surprisingly, breezes instead 
(jKasting fc Pollacklll983HYelle.2004l . Relaxation meth- 
ods are also suitable for our problem because the wind 
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profile is expected to be smooth, with no osciUatory be- 
havior, and so pre-defining a radial grid for the solution 
is not especially problematic. Nevertheless, care needs 
to be taken in implementing the method; we describe 
as follows our procedure, developed after considerable 
experimentation. Sections 12.2.11 — [2.2.41 describe how we 
compute the wind profile from some base depth in the at- 
mosphere up to the sonic point; section [2.2.51 takes this 
solution and extends it out to the Ro che lobe. 

F or alternate num erical methods, see lTian et al.l ()2005[ ) 
and lGarcfa Mufiozi (^07). 

2.2.1. Finite Difference Equations: From the Base to the 
Sonic Pomt 

From a location r^-m = Rp in the upper atmosphere 
of the planet to the sonic point Ts of the wind, we 
use the Numerica l Recipes relaxation routine solvede 
(jPress et al.lll992l ) to solve the finite difference versions 
of©, (ED, ©, ®, and 0: 



E 



= A,p 
= A,p- 



E- 



2j 



E: 



3j ■ 



:A,T- 



^Ar 

dr ' 

p[- + --r 

\r V dr 
dv . 
dr 



— "fkT/ II 



A,r = 



(9) 



27fcr/(H - (7 - 



A,r = 



(10) 
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\pv k 
T df+ 



Tdp 
p dr 



(1 + /+) dr 



A,r 







Eij = AoT Ao7 

dr 



= Af, 



.oAj-r = 



(11) 



(12) 



E,,^A,U^^Ap 



pv 



uve 



Up 



(1 - f+)P 



A,r = 



(13) 

where Q = T + A, and A^a; = Xj — Xj^i at the jth radial 
grid point. In evaluating the individual variables that 
make up the derivatives, we average across adjacent grid 
points; e.g., p = (pj + o,-i)/ 2. This is the same choice 
adopted bv [Press et all (|1992l ). 

We introduce an extra dependent variable z = rg — rmin 
because we do not know a priori the location of the sonic 
point Tg. Thus 



E^ 



6j 



A,z 







(14) 



In other words, we solve for z just like we do any other 
dependent variable, and its solution tells us the radial 
location of the sonic point (rg). 



Our six dependent variables are p, u, T, /+, r, and 
z, all non-dimensionalized for ease of calculation us- 
ing the scales po = 10~^^ g/cm^, Tq ^ lO"' K, vq = 
(fcTo/TOH)^/^ = 9 km/s, and zq = 10^" cm. We solve 
for these variables on a radial grid of 1000 points that 
has more points concentrated near rmin (where deriva- 
tives are large) and near Tg. The convergence parameter 
"conv" that solvede uses is set to 10~^°. The scale pa- 
rameters "scalv" used to calculate the convergence pa- 
rameter are p = 100, i; = 2, T = 1, r = 100, /+ = 1, 
and z = 3 in our non-dimensionalized units. We take 
our independent, radius- like variable to be q such that 
r — faiin + qz; q runs from to 1. 

The finite difference equations are solved by a multidi- 
mensional Newton's method, which requires that we eval- 
uate partial derivatives of the Eij's with respect to the 
dependent variables. We evaluate these partial deriva- 
tives numerically, by introducing small finite changes in 
each of the dependent variables and calculating the ap- 
propriate differences. 

2.2.2. Boundary Conditions 

We need six boundary conditions (BC) to solve Equa- 
tions ©-(mi). Two boundary conditions are provided by 
the requirement that the wind pass through the critical 
point where 



(BCl) 



and 



2-fkT GMp (7-l)Qr 3GAf* 



P 



pv 



(15) 



= (BC2) 



(16) 



in order to avoid infinite derivatives (see Equatio n [TO 



These are the critical point cond itions of the iParker 
(1958) transonic wind (see, e.g., iLamers fc Cassinelh 
119991) . 

We choose our remaining four boundary conditions 
as follows. At the base of the fiow, we set p(rinin) = 
4 X 10-"gcni-3 (BC3); f+{r,^in) = 10"^ (BC4); and 
T — 1000 K (BC5), approximately the eff ective tempera- 
ture of a hot Jupiter at a = 0.05 AU (e.g.. lBurrows etld] 
I2OOI . For our final boundary condition, we enforce the 
condition that r(rs) equals the optical depth between 
the sonic point and the Roche lobe — see §2.2.51 For 
our standard model, this optical depth turns out to be 
r(rg) = 0.0023 (BC6). 

In Appendix [Aj we demonstrate that our solution is 
insensitive to these particular choices of numbers for 
BC3 through BC6. We show there that the solution 
hardly changes as long as p(r-n-im) is large enough that 
T(rinin) > 1; /+(r-min) < 1; ^(rmin) <C lO'^ K; and 
r(rs) <C 1. We also describe how the mass loss rate 
changes by less than a factor of 2 for 10% variations in 
T-m\\\ about our fiducial radius i?p = 10^" cm (see the 
order-of-magnitude discussion at the beginning of ij2l for 
why Tijiin is only uncertain by about 10%). The insensi- 
tivity to boundary conditions helps to justify our neglect 
of the II2/II dissociation front, which is located at greater 
depth than the II/II+ ionization front — the latter we do 
resolve, near rmin. 
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2.2.3. Sonic Point Limit 

Because the exact expression for dv/dr in Equation 
PI}|) — and, by extension. Equations ^ and pT|) — is dif- 
ficult to evaluate accurately near the sonic point (both 
numerator and denominator of dv/dr vanish there), we 
have derived an analytic form for it that is strictly valid 
only at the sonic point: 



dv 
dr 
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(17) 



where Qi and Q2 are defined by dQ/dr = Qi{dv/dr) + 
Q2- This expression is derived by applying L'Hopital's 
rule to dv/dr in Equation 1 10( and simplifying the result 
using boundary conditions (fT5|) and 
In our code, we evaluate dv/dr as 



dv dv 

~r~ — exact ~r~ 
dr dr 



where 



— —erf 



P 



(1-F_0^ 
dr 



^ lkT\ 



(18) 



pv^ 



erf is the error function, p = 100 is a parameter that 
determines the width of the transition between the two 
right-hand terms, and (dw/(ir)oxact is given by Equation 
pU|) . Far from the sonic point, {dv/dr)r^ is not accurate, 
so where inexact is within machine precision of 1 , we revert 
to using (dw/dr)cxact only. 

2.2.4. Solving Successtvely More Complicated Problems 

Relaxation codes require good initial guesses to con- 
verge, so we build our final solution by solving succes- 
sively more complicated problems. The solution of a 
given problem furnishes the initial guess for the next 
problem. 

First we use our relaxation code to find solutions for 
the separate problems of an isothermal wind without 
photoionization, and an isothermal hydrostatic atmo- 
sphere with photoionization, both neglecting tidal grav- 
ity. The combination of these solutions provides the ini- 
tial guess for an isothermal wind with photoionization. 
Then we remove the restriction that the wind be isother- 
mal. The photoionization heating term, followed by cool- 
ing terms and finally the tidal gravity term, are added 
one by one. Sometimes the addition of even a single term 
requires iteration: the term must be added in diluted 
form first and gradually strengthened to full amplitude. 

2.2.5. From the Sonic Point to the Roche Lobe 

From the sonic point outward, we use the Numerical 
Recipes routine od eint with a Bulirsch-Stoer integrator 
(jPress et al.l 11992*) to solve our original ordinary differ- 
ential equations ('i J2.ip . We set the routine's convergence 
parameter "eps" to 10^^'^. The solution is extended out 



to the planet's Roche lobe radius i?Rochc- This extension 
is a straightforward initial value problem, with initial 
conditions at the sonic point provided by our relaxation 
solution. For the first few radial steps in the integration, 
we use ()17p for dv/dr to avoid numerical singularities. 

Finally, this extended solution feeds back into our re- 
laxation code through T(rs). We iterate a few times be- 
tween the relaxation code and the Bulirsch-Stoer integra- 
tor until T(rs) self-consistently reflects the optical depth 
between the sonic point and the Roche lobe radius. The 
assumption here is that the optical depth beyond the 
Roche lobe is negligible; this assumption is valid because 
jirs) ■€i 1 for solutions of interest to us (see Appendix 

2.3. Results 

In i i2.3.H we present the results for our standard 
model, appropriate for a hot Jupiter orbiting a Sun-like 
star on the main sequence. In ii2.3.2[ we present a sample 
high flux case for which F\jy is increased a thousandfold 
over its standard value, as would be the case for a hot 
Jupiter orbiting an active pre-main-sequence star. In 
§2.3.31 we describe how the maximum mass loss rate M 
varies when -Fuv ranges over four decades. 

2.3.1. Standard Model: Hot Jupiter Orbiting Main-Sequence 

Star 

Figures [TH5] display the results for our standard model. 
This numerical solution verifies many of the order-of- 
magnitude estimates made at the beginning of [J2l The 
stellar UV flux drives a transonic wind with temperature 
^10"* K and velocity ^^10 km/s. The hydrogen density 
where optical depth unity to Lyman continuum photons 
is reached is ^lO^cm"'^. Because our solution pertains 
to the substellar ray connecting the planet to the star, it 
yields the maximum mass flux pv; escape is most aided 
by tidal gravity along this ray, and the substellar point 
receives the greatest UV flux. Applying this solution over 
the entire surface of the planet yields an upper bound on 
the mass loss rate of M = 3.3 x 10^° g/s. In gS^HSSl we 
estimate the factors by which the actual mass loss rate 
is reduced. 

Figure [3] (bottom panel) displays the relative contri- 
butions to ionization balance as a function of altitude. 
Above the r = 1 surface, gas advection, not radiative re- 
combination, balances UV photoionization. A hydrogen 
atom, once photoionized, does not have time to recom- 
bine before it is swept outward with the wind. As a gas 
parcel travels outward, its electron density decreases and 
the recombination time l/(n+Q!rcc) becomes ever longer; 
more and more of its atoms are ionized, and the ioniza- 
tion fraction increases with altitude. At the sonic point, 
about 20% of the hydrogen remains neutral. This situ- 
ation differs from static H II regions in which photoion- 
ization is balanced by radiative recombination and the 
transition between ionized and neutral gas is sharp. 

Figure [3] (top panel) also shows that the temperature 
of the gas is largely controlled by heating from photoion- 
izations and by cooling from gas expansion (PdV work), 
with a small contribution from cooling by Lya radiation 
at the wind base. In the upper portions of the wind, 
cooling by PdV work lowers the gas temperature from 
its peak of 10000 K to about 3000 K (Figure [l]). Since 
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Fig. 1. — Standard planetary wind model, which uses param- 
eters inspired by the hot Jupiter HD 209458b. A UV flux of 
450ergcm~^ s~^ incident on a planet with mass 0.7Mj and radius 
Rp = 1.4i?j located 0.05 AU from a IMq star drives a transonic 
wind. The flow is calculated along the line joining the planet to 
the star. Density p (top), wind speed v (middle), and temperature 
T (bottom) are presented as functions of altitude. On each panel, 
the sonic point of the wind is marked with an x, and the t = 1 
surface to photoionization is marked with a diamond. The planet 
loses mass at a maximum rate of M = Aizr^pv = 3.3 X 10^*^ g/s. 

these temperatures exceed the '^2000 K required for H2 
to dissociate, our neglect of H2, and by extension the ra- 
diative coolant that is known to be important wher e 
H2 exists in abundance ()Yellell200^ iGarcia Mufiozll200l . 
is self-consistent, at least above the r = 1 surface to pho- 
toionization of atomic H. In other words, most of our 
modeled region is self-consistently devoid of molecular 
hydrogen. 

Integrated over the entire radial extent of our model 
from Tniin to i?Rochc7 & mcasurc of the heating rate per 
solid angle (measured from the planet) from photoioniza- 
tions is Lphoto = / Fr^dr = 2.3 x 10^^ ergs~^ sr~\ about 
equal in magnitude to the analogously defined cooling 
rate due to PdV work, Lpdv — —1-9 x 10^^ ergs~^ sr~^ 
In comparison, the height-integrated cooling rate from 
Lya radiation is only Llyq = —2.9 x 10^^ ergs~^ sr~^. 
(Note that iphoto+ipdv +-^Lya 7^ because the internal 
energy of the gas changes along the flow; see Equation[3|) . 
The fact that as much as 83% of photoionization heat- 
ing goes into PdV work implies that mass loss is largely 
"energy-limited"— see ij2. 3.31 and gSH 

We have used our solution to estimate, ex post facto, 
the contributions to internal energy changes from recom- 
bination radiation, free-free emission, coUisional ioniza- 



r/R„ 



Fig. 2. — Planetary wind model for the same standard parame- 
ters presented in Figure[T] Number densities of neutral plus ionized 
hydrogen n, neutral hydrogen no, and ionized hydrogen n+ (top); 
ionization fraction = n^/n (middle); and optical depth t to 
photoionization (bottom) are presented as functions of altitude. 
On each panel, the sonic point of the wind is marked with an x, 
and the r = 1 surface is marked with a diamond. As much as 
--^20% of the hydrogen remains neutral at high altitude. 

tions, and conduction. As demonstrated in Figure [3l 
these are not significant. 

2.3.2. High UV Flux Case: Hot Jupiter Orbiting a T Tauri 

Star 

Figures HHS] display results for a much larger UV flux of 
FuY = 5 X 10^ erg cm~^ s~^, characteristic of the radia- 
tion field e xperienced by a hot Jup iter orbiting a T Tauri 
star (e.g.. ' HoUenbach et all 120001 ). A gain, the UV flux 
drives a transonic wind with temperature ~10'* K and 
velocity a few x 10 km/s. The density of the wind is sub- 
stantially greater, however, than in the standard model. 
Applying our solution over the entire surface of the planet 
yields a maximum mass loss rate of M = 6.4 x 10^^ g/s. 

In contrast to the standard model, here radiative re- 
combination, not gas advection, balances photoioniza- 
tion above the r = 1 surface. The transition from a 
neutral to a nearly fully ionized flow is sharp. Also in 
contrast to the standard model, Lya cooling plays a dom- 
inant role in balancing photoionization heating near the 
base of the wind. The global energy budget now divides 
as follows: Lphoto : Lpdv ■ Li^ya = 1 : —0.29 : —0.67 
(again, iphoto + Lpdv + ^Lya ^ because the internal 
energy of the gas changes along the flow) . Thus the high 
flux case is strongly radiatively cooled. These properties 
of the high flux case are similar to those of classic H II 
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photoionization heating 



Lyoc cooling cooling by 
PdV work 




photoionization 

gas advection 




Fig. 3. — Contributions to energy equilibrium (top) and ion- 
ization equilibrium (bottom) for the standard transonic wind dis- 
played in Figures [T] and [2] In each panel, the sonic point of the 
wind is marked with an x, and the r = 1 surface to photoionization 
is marked with a diamond. Energy balance is dominated by pho- 
toionization heating and PdV work done by expanding gas, with 
some contribution from Lyo cooling near the base of the wind. 
Internal energy changes due to recombination radiation (r.r.), col- 
lisional ionization (c.i.), free-free radiation (f.f.), and conduction 
(either cooling c.c. or heating c.h.) are negligible (see Appendix A 
for the formulae used to evaluate these extra contributions). UV 
photoionization is balanced by gas advection above the r = 1 sur- 
face. 
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4. — Same as Figure [T] except that -Fuv = 5 X 
~^ s~^ . Such a flux is characteristic of that incident on a 



Fig. 
10^ erg cm 

hot Jupiter orbiting an active pre-main-sequence star. The planet 
loses mass at a maximum rate of M = inr^ pv = 6.4 X 10^^ g/s. 



variable tidal gravity. In H'SA[ we consider how the plan- 
etary outflow interacts with the ambient stellar wind — 
and how such interaction can strongly suppress dayside 
winds while energizing nightside winds. Finally, in §3.5( 
we assess whether stellar radiation pressure can drive sig- 
nificant planetary outflows. 



regions (Stromgren spheres). Mass loss under these con- 
ditions is "radiation/recombination- limited" — see §2.3.31 
and g321 

2.3.3. Mass Loss Rate vs. UV Flux 

We calculate the maximum (i.e., spherically symmet- 
ric) mass loss rate as a function of incident UV flux 
and display the result in Figure [71 For Fuv less than 
~10 erg cm" s"\ M oc F^-^. For larger UV fluxes, the 

mass loss rate increases more slowly as Af oc F^y. We 
discuss the origin of this difference in 

3. DISCUSSION 

In we check the validity of our assumption that 
the flow is hydrodynamic. In §3.2) we explain how there 
exist two modes for the wind, and derive analytic expres- 
sions for the mass loss rate that reproduce fairly well our 
numerical results. In i i3.3l the assumption of spherical 
symmetry is relaxed; we estimate the factors by which 
mass loss rates are reduced because of variable UV irra- 
diance across the surface of the planet, and because of 



3.1. Hydrodynamic vs. Jeans Escape 

We have modeled the mass outflow from a hot Jupiter 
as a hydrodynamic wind rather than as local Jeans es- 
cape. For this fluid description to be accurate, the 
gas must remain coUisional until it is above t he sonic 
point of the flow (jChamberlain fc Hunteiilll987l p. 377). 
In other words, the exobase — the height at which the 
scale length H = p{dp/dr)~'-^ equals the mean free 
path to collisions Amfp — l/(ncrpp), where Cpp ~ 
10~^'^(T/10^K)~^ cm"^ is the Coul omb cross sec tion for 
protons scattering off protons (e.g.. lSpitzeilll978h — must 
lie above the sonic point. This requirement is easily 
satisfied by our models. For our standard model with 
FiJV = 450ergcm~^ s~-^, if/Amfp ~ 10* at the sonic 
point. For Fuv = 5 x 10^ ergcm"^ s~\ H/Xmfp ~ 10^ at 
the sonic point. Therefore the gas behaves as a collisional 
fluid, and mass loss rates b ased on Jeans escape cri- 
teria (|Lecavelier des Etangs e t al. 2004; Ycllc 200|) are 
not appropriate. We agree with ITian et al.l (|2005l ) and 
iGarcfa Muhod (|2007f ) that mass loss from hot Jupiters 
is hydrodynamic, not ballistic. Figure [8] summarizes the 
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Fig. 5. — Same as Figure [2] except for the high UV flux case. 
By contrast to the standard model, the wind here is more fully 
ionized, and transitions from neutral to ionized more sharply. 

scales in our wind solutions. 

Note that for our standard model of the substellar 
streamline, the sonic point lies inside the Roche lobe, 
and therefore our account of tidal gravity in ^ is ade- 
quate. There may be other streamli nes where th e sonic 
point is outside the Roche lobe. iJaritz et al.l ()2005f ) 
suggest that for such streamlines, mass loss occurs via 
"geometrical blow -off" and t he Jeans escape criteria 
of Lccavelicr dcs E tangs et al.l ([2004) should be applied. 
We disagree; even if the sonic point lies outside the Roche 
lobe, gas pressure gradients inside the Roche lobe will 
drive an outflow, and wherever the gas remains coUi- 
sional, the flow must be solved using the equations of 
hydrodynamics . 

3.2. Energy- Limited vs. 
Radiation/Recombination- Limited 

The two regimes for mass loss that our numerical so- 
lution uncovered — M oc F^y at low Fuv and M oc F^y 
at high Fuv (Figure [7]) — can be understood simply. 

At low Fuy, the flow is largely "energy- limited." Most 
of the energy deposited by photoionizations as heat — 
i.e., '^e-TrFuv-Rp — goes into PdV work, with little loss to 
radiation and internal energy changes (the relative con- 
tributions to the energy budget are given in ij2.3.f p . The 
PdV work lifts material out of the planet's gravitational 
potential well: measured per unit mass, the work done is 

PAV PRl pgHRl GMp 
pR^H ^ pR^H ^ pR^H ^ Rp ' 




Fig. 6. — Same as Figure [S] but for the high UV flux case. By 
contrast to the standard model, Lya cooling balances photoioniza- 
tion heating at the wind base, and radiative recombinations balance 
photoionizations everywhere. The wind here behaves more nearly 
like an standard H II region. 



Then the energy-limited mass loss rate is given by 

V0.3/ VIO °cmy V ^^P / 

(450ergcm-2s-i) ^^^^ 

close to the result found numerically at low -Fuv (the 
factor of 5 difference in normalization between Equation 
I19l and the curve shown in Figure [7] arises mostly because 
the latter takes the substellar UV flux and applies it over 
all 4tt steradians, whereas the former averages the UV 
flux over the surface of the planet — hence the factor of tt 
in equatioii ll9D . E nergy - limited outflows were also found 
by I Watson et al.l (|l98ll ) who studied mass loss in the 
highly conductive atmospheres of the terrestrial planets.^ 
At high -Fuv I the flow is "radiation/recombination- 
limited." As quantified in i )2.3.2i the input UV power 
is largely lost to cooling radiation. Radiative losses ther- 
mostat the gas temperature to T ~ 10"* K. Under the ap- 
proximation that the wind is isothermal, M ~ AnpsCsr^ 

^ Note, however, that IWatson et al.l l|19811 ) reserve the phrase 
"energy-limited" for use in another context. Neverthel ess their 
Equation (2) is essentially the same as our Equation JToJl, the 
"energy-limited" mass loss rate in the sense that we use the phrase. 
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Fig. 7. — Maximum mass loss rates as a function of UV flux for a 
0.7A/j planet with radius 1.4i?j located 0.05 AU from a 1A/q star. 
We calculate the mass loss rate along the substellar ray connecting 
the planet to the star, along which escape is most easily achieved, 
and use this solution over the entire surface of the planet to cal- 
culate the maximum M = Airr'^pv. More realistic estimates that 
average the incident UV flux over the planet ary s urface yield mass 
loss rates that are lower by factors of 3-4 f t|3.3l l. For Fuv ^ 10^ 
erg/cm^/s, M oc F^y- For Fuy > 10^ erg/cm^/s, the mass-loss 
rate increases more slowly as M oc F^y. These two regimes corre- 
spond to "energy-limi ted" and "radiation/recombination- limited" 
flows, as explained in i|3.2l 



at the sonic point, where Cg = [kT/{mii/2)]^^'^ is the 
isothermal sound speed (the factor of 2 accounts for 
the fact that the hydrogen is nearly completely ionized) , 
= G'Mp/(2Cs), and ps is the sonic point density. -'^^ Be- 
tween the r = 1 surface and the sonic point, the density 
structure is nearly hydrostatic, so that 



Pbaso exp 



GM„ fR 



R r2 



1 



In the high flux case, the density pbaso at the t — I 
photoionization base is n_|_,basc"^H- That density is de- 
termined by ionization equilibrium, which involves a bal- 
ance between photoionizations and radiative recombina- 
tions (Figure [6|): 



UV 



■CiyoJT-O.basc ^ '^4^,bascQ^i-cc , 



the 



neglecting the order unity flux attenuation at 
T = 1 surface. The base neutral density np^baso ~ 
l/{(JvgHbasc) ^ mug / {2auokT) . This neutral density is 
fairly insensitive to F\jy (see beginning of [J2] and com- 
pare Figures O and and so we conclude that the 
radiation/recombination-limited mass loss rate is given 

by 



M„-lim ~ 4 X 10 
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UV 



5 X 10^ erg cm ^ s i 



1/2 



(20) 



similar to the answer found numerically at high i<uv- 

See, e.g.. lLamers &: Cassinelii|[T999l . 
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Fig. 8. — Summary of important scales in our wind solutions. 
For a planet with properties inspired by HD 209458b, the t = 1 
surface to which UV photons penetrate is located where the local 
pressure is measured in nanobars, at about 1.1 times the radius 
of the 1 bar surface. Below this photoionization base, the temper- 
ature is ^lO"^ K, close to the effective temperature, as befits an 
atmosphere in radiative equilibrium with the incident stellar opti- 
cal radiation balancing the planet's thermal emission. Molecular 
chemistry (which we do not model) is important in this regime. 
Our model is valid above the r = 1 surface, where the temperature 
is thermostated to ~10* K by Lyman-a cooling. We demonstrate 
in Appendix |A] that our standard wind solution is insensitive to 
conditions below the t = 1 surface, justifying our simple model 
which treats only atomic and ionized hydrogen. The sonic point of 
our wind solutions lies at 2-4 planetary radii, with smaller sonic 
point distances corresponding to higher UV fluxes. The Roche lobe 
radius of the planet at '~4.5iJp is close enough to the sonic point 
radius that tidal gravity is significant. Below the exobase, where 
the mean free path to collisions equals the scale height of the at- 
mosphere, the gas behaves as a fluid. Because the exobase is well 
outside the sonic point of our winds, mass loss from hot Jupiters 
takes the form of a hydrodynamic wind rather than Jeans escape. 

3.3. Spherical Asymmetry: Day/Night and Tidal 
Gravity 

The mass loss rates given in all our plots are upper 
limits because they take our ID solution for the sub- 
stellar streamline and apply it over 47r steradians. The 
mass flux is maximized for the substellar streamline be- 
cause the substellar point receives the maximal UV flux, 
and because the tidal gravity term weakens the planet's 
gravity most along the line joining the planet to the star. 
We now discuss the extent to which the actual mass loss 
rate is reduced because of day/night differences in the 
received UV flux, and because of the directional depen- 
dence of tidal gravity. For simplicity we discuss these 
effects separately, as if they could be isolated from one 
another. 

Ignoring for the moment differences between day/night 
external boundary conditions — these are actually signifl- 
cant and dealt with in H3.4I — day/night differences would 
be erased in the extreme case that horizontal winds re- 
distribute photoionized plasma from the dayside of the 
planet to the nightside, on a timescale shorter than the 
wind's radial advection time of a few hours. In this case, 
the "nightside wind" would blow just as strongly as the 
"dayside wind." In reality, if the dayside wind blows 
freely (see t j3.4l for important reasons to believe that it 
may not), the timescale for horizontal advection is at 
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least as long as that for radial advection, since the dis- 
tances travelled in both cases are several Rp, the radial 
speed is supersonic, and the horizontal speed is at most 
sonic. So if the dayside wind blows, the nightside wind 
is likely muted. The inverse is also true — see §3.41 

In the opposite limit of no horizontal redistribution, 
the mass loss rates we have calculated would be reduced 
by a factor of (1/2) J^^^ sm9{cos9)''de = l/[2(7 + 1)], 
where 9 is the angle measured from the substellar ray 
{9 = points along the substellar ray, while 9 = tt/2 
defines the terminator dividing day from night), and 
the (cos^)''' factor accounts for how the planetary mass 
flux scales with incident stellar flux (which itself scales 
as cos 9). For energy-limited flows, 7 = 0.9, and for 
radiation/recombination- limited flows, 7 = 0.6 (see Fig- 
ure [7]). The reduction factor equals 1/3.8 = 0.26 and 
1/3.2 = 0.31 in the two respective cases. If there is hor- 
izontal redistribution, we do not expect the reduction 
factors to change appreciably from these values. To first 
order, redistribution of plasma simply redistributes the 
wind over the planetary surface. Changes in the total 
mass loss rate are expected to be of second order. For 
example, if the wind is strictly energy- limited fS )3.2p . the 
mass loss rate depends only on the amount of UV radia- 
tion intercepted by the planet, and is independent of the 
degree of redistribution. 

To get a sense of how much the mass loss rate is re- 
duced because tidal gravity does not point parallel to all 
streamlines emanating from the planet, we eliminate the 
tidal term from the force equation ([2]) and re-solve the 
fluid equations. We obtain a mass loss rate (multiplying 
by 47r steradians) that is 0.79 x our standard model re- 
sult. This can be considered a maximum reduction factor 
for the directional dependence of tidal gravity insofar as 
its effect can be isolated. 

Though tidal gravity does not change the calculated 
mass loss rate appreciably, it is nonetheless important 
for two reasons. First, tidal gravity alters the velocity 
structure of the wind, allowing it to accelerate to larger 
speeds (Figure [9]). Second, tidal gravity moves the loca- 
tion of the sonic point to lower altitude, allowing it to 
remain within the Roche lobe along the substellar ray 
and helping validate our ID treatment. 

To summarize this subsection, if the dayside transonic 
wind blows freely, our crude accounting for the direc- 
tional dependences of UV irradiance and of tidal gravity 
suggest that actual mass loss rates are lower than our 
plotted values by factors of ~4. We now turn to the 
question of whether the dayside wind can actually blow, 
when faced with the considerable pressure of the host 
star wind. 

3.4. Colliding Winds and Breezes: Trading the Dayside 
Wind for the Nightside Wind 

The planetary wind does not exist in vacuum. The 
dayside wind blows into the plasma streaming from the 
star: the stellar wind. Stellar and planetary winds col- 
lide and mix in a standing bow shock surrounding the 
planet. The situation is a nalogous to the colliding winds 
of m assive stellar binaries (|Luo et al.lll990l : IStevens et al.l 
|199 2[) . Figure \TU\ supplies a cartoon illustration. How- 
ever, as stressed in the figure caption and in the discus- 
sion below, numerous assumptions underlie this cartoon, 
and reality is likely to look substantially more compli- 
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Fig. 9. — Wind velocity as a function of altitude for our standard 
model with tidal gravity (solid curve) and without (dashed curve). 
Tidal gravity alters the velocity profile of the wind, allowing the 
flow to attain higher speeds. It also moves the wind's sonic point 
to lower altitude. 
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Fig. 10. — Cartoon of the standing bow shock (heavy curve) 
where planetary and stellar winds collide and mix. For this figure 
only, in order to compute the location of the bow shock, we assume 
for simplicity that the winds are spherically symmetric, have zero 
thermal and magnetic pressure, and have constant velocity (we 
neglect Coriolis, centrifugal, and gravitational forces — wc have not 
directly used one of our wind models). For the stellar wind, we 
assume a mass loss rate of 5 X 10^^ gs~^ and a velocity of 400 
km/s. For the planetary wind, we take 3 X lO^'' gs"! and 20 km/s, 
values informed by our wind models. The bow shock is given by 
the condition that the bulk mom entum fluxes ("r am" pressures) 
normal to the shock balance (e.g.. lLuo et al.lll990l) . Such a shock 
does not influence the planetary wind upstream if it is located 
outside the sonic radius, as is drawn. In reality, this condition may 
not be satisfied because of the overwhelming total pressure of the 
stellar wind. 

cated. 

At the stand-off shock, the total pressures of the two 
flows balance. The total pressures include contributions 
from "ram" (pw^), thermal, and magnetic pressure, but 
should include only components that are normal to the 
shock front (the shocks are oblique). The total pres- 
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sure of the planetary wind decreases monotonically with 
distance from the planet. If the total stellar wind pres- 
sure P*,tot near the planet is less than the total plane- 
tary wind pressure Ptot(''s) at the planetary wind's sonic 
point, then the stand-off shock will be located down- 
stream of the planet's sonic point. There it cannot in- 
fluence the planetary flow upstream, for the same reason 
that one cannot shout upstream in a supersonic flow and 
be heard. In this case the transonic wind solution that 
we have obtained for r < Vs — and in particular our com- 
puted mass flux — would remain unchanged. 

What are P*,tot and Ptot('^s) for a hot Jupiter orbiting a 
main-sequence, solar type star? For the former quantity 
we are guided by observations of the solar wind. At a = 
0.05 AU = lOi?0, th e solar wind h as accelerat ed to nearly 
its m aximum speed ([McKenzie et al. 1997: Marsch et al.l 
|2003() . Approximating the stellar wind velocity as con- 
stant from 0.05 to 1 AU, we scale characteristic space- 
craft measurements (SOHO^^; Ulysses^'^) of the solar 
wind density and velocity at 1 AU to estimate the cor- 
responding density and velocity at 0.05 AU: n, ~ 
6(1/0.05)^ protons cm~^ and ii* ~ 400kms~^. We take 
the local proton temperature to be T* ~ 10^ K (the elec- 
tron temperature is several times lower; iMcKenzie et al.l 
|1997() . and the magnetic field stren gth to be ~ 0.01 
G at 10 stellar radii ( Banaszkicwic z et al.lll99^ . Then 
-P*,tot ~ f^*WHW* +n^,kT^, + B'^/{8tt) ~ 10 picobars, with 
the ram and magnetic pressures dominating. This is an 
overestimate insofar as we are not taking the compo- 
nents of the ram and magnetic pressures that are rele- 
vant perpendicular to the shock front. Compare this to 
-Ptot(?'s) ^ 3 pbars, as calculated for our standard model 
of the substellar streamline (see Figure fTTj). neglecting 
the unknown magnetization of the planet. Given the 
uncertainties, the most we can say is that P*.tot and 
-Ptot(?'s) are comparable. While a transonic wind solution 
along certain (not necessarily substellar) streamlines may 
yet be possible, it is also possible — indeed even likely, 
given the fact that the solar wind can gust to values of 
and several times higher than the ones we have 
used — that the stellar wind squashes the planetary day- 
side outflow down to a subsonic breeze, or even stops 
dayside photoevaporative mass loss completely. 

A breeze would still blow a bubble in the stellar wind, 
like that drawn in Figure [TUl only smaller. But because 
a breeze is subsonic, it would not traverse a shock at 
the bubble boundary. Unlike the case for the super- 
sonic wind, a breeze is in causal contact with the flow 
at the boundary; i.e., conditions at the boundary influ- 
ence, via sound waves, conditions in the interior. The 

Though magnetic disturbances can propagate upstream if the 
stand-ofT shock occurs inside t he planetary wind's Al fvcnic and fast 
magnetosonic points (see, e.e.. lWeber fc DavisillQGTl for a theory of 
magnetized winds). 

http://umtof.umd.edu/pm/ 

http:/ /swoops. lanl.gov/recentvu. html 

It has been speculated that the magnetic fields of hot Jupiters 
are weaker than that of Jupiter, since the rotation of a hot Jupiter 
is tidally slowed to nearly the orbital period of three days (exact 
synchronization is not possible if the planet lacks a permanent 
quad rupole moment; see, e.g., the textbook by [Murray &: Dermotl 
I200C1I| . whereas Jupiter's rotation period is ten hours. iBlackett 
11194 7) speculates that planetary magnetic moments scale linearly 
with rotation frequency; if so, the surface field on a hot Jupiter at 
r = Rp would be G. If the field falls as a dipole to r = rs ^ 3-Rp, 
it would add of order 10 pbar to Ptot(''s)- 



radial velocity would decrease smoothly to zero at the 
bubble boundary, giving rise to large pressure gradients 
that would drive flows parallel to the boundary (by anal- 
ogy to subsonic flow around a blunt obstacle such as a 
hard sphere). The flow inside the bubble would be non- 
radial and multidimensional. Unfortunately, we cannot 
easily capture such behavior with our ID model; if we 
tried to model a breeze by imposing a boundary con- 
dition of zero radial velocity at flnite distance from the 
planet, while insisting simultaneously on a non-zero mass 
loss rate, our simple continuity relation ([1]) would yield 
inflnite density. 

Nevertheless, we can get a sense of how much the day- 
side mass loss rate might be reduced in the presence of 
external pressure, by re-running our relaxation code with 
boundary condition replaced by [u^ = f3jkT/ fi]r^, 
where /3 < 1 is our breeze parameter. Note that rg in this 
case is no longer the sonic point but is instead the loca- 
tion where the flow stops accelerating and starts deceler- 
ating. Resulting breeze solutions are shown in Figure [TTl 
with the corresponding mass loss rates plotted in Figure 
fT2l While the total pressure Ptot for the (3=1 wind so- 
lution decreases to zero at large distance, Ptot for a /3 < 1 
breeze asymptotes to a finite value. As /3 decreases, the 
breezes blow more slowly, approaching the v = limiting 
solution for a hydrostatic atmosphere. Once (3 < 10~^ 
{v < O.ly/jkT/^), the Ptot profiles hardly change. We 
take as a reference pressure Ptot('"s); the location is 
a sensible one to examine since for the most part it de- 
creases as P decreases, mimicking the shrinking of the 
planetary bubble with stronger stellar winds. Over the 
entire family of breezes, the reference pressure Ptot(?'s) 
achieves a maximum value of ~40 pbar for our standard 
parameters and boundary conditions. (Recall that an 
isothermal hydrostatic atmosphere has a finite pressure 
at infinity.) 

To the extent that our ID solution lends insight into 
the 3D breeze problem, we conclude the following: if the 
stellar wind pressure P*,tot ^ 3 pbars, the planet emits 
a full-fledged transonic wind; if P*,tot lies between '^3 
and ~40 pbars, then the planet emits a dayside breeze 
with a flnite mass loss rate; but if P*^tot ^ 40 pbars, 
then the planet's dayside atmosphere is forced to be ra- 
dially hydrostatic, with the stellar wind penetrating to 
the depth where pressure balance is achieved. Figure 
1121 informs us that tiny increases in Ptot(fs) above 35 
pbars correspond to enormous reductions in M. At the 
same, if Ptot(T's) is below 35 pbar, the mass loss rate is 
^10^" g s~^, to within a factor of about 3. Thus, the stel- 
lar wind can act essentially as an on-off switch — at times 
permitting the planet to lose mass from its dayside at 
near maximum rates of order 10^°gs~^, and at others 
shutting down dayside mass loss completely. In the case 
of hot Jupiters orbiting main-sequence stars, given how 
the stellar outflow pressure is of order 10 pbar and how 
it may increase dramatically during violent coronal mass 
ejections, we expect this switch to flip back and forth 
with stellar activity. We cannot rule out the possibility 
that the dayside switch might even be off the majority 
of the time. 

Even when the dayside switch is off, however, we do not 
expect mass loss to shut down completely. As the dayside 
outflow weakens, the nightside wind should strengthen 
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Fig. 11. — Velocity and total pressure (Ptot = p'o^ + pro- 
files for our unique transonic wind (/3 = 1) and a family of breezes 
(/3 < 1 as labelled), for standard model parameters. The x sym- 
bol marks r = rs, which for the wind represents the sonic point, 
and for breezes represents the point of maximum velocity. At in- 
finity, breezes remain pressurized while the wind does not. In the 
middle panel, profiles not labelled include fi = 10~^, 10~^, and 
10~*. Note how these profiles hardly change as /3 decreases — the 
solution is approaching that of a d = hydrostatic atmosphere. 
The bottom panel shows that for the wind, thermal pressure dom- 
inates at depth (dotted line), while pv^ "ram" pressure dominates 
at altitude (dashed line). For the /3 = 0.01 breeze shown, thermal 
pressure dominates everywhere. For the most part, shrinks as 
ever slower breezes (smaller jS) are considered. The exception oc- 
curs for /3 < 10~^ (middle panel), and is an artifact of our choice 
to compute T(rs) out to the fixed Roche lobe radius -RRochci ''"(^s) 
increases with decreasing /3 < 10^^. 

in proportion. That is because the UV energy that the 
planet absorbs must be spent one way or another; if it 
is prevented from doing PdV work against gravity on 
the dayside, it wiU instead power horizontal winds which 
carry heat to the nightside."'^^ Shielded from the stellar 
wind, the nightside is free to emit its own wind. Thus, 
the total mass loss rate might remain roughly constant at 
'^lO^^gs"^, even when the stellar wind gusts (assuming 
that no other energy sink becomes active). 
In concluding that nightside winds will blow when day- 

IKoskinen et al.l COOyal lbt) construct a thermospheric circula- 
tion model for hot Jupiters similar to HD 209458b but located 
further than 0.16 AU from their host stars. They argue that at 
these distances, atmospheric temperatures remain low enough that 
molecular hydrogen never dissociates, Hg" dominates cooling, and 
the atmosphere remains in vertical hydrostatic equilibrium. Sim- 
ilar circulation models could be made for the thermospheres of 
atomic hydrogen of closer-in planets, in the case where the dayside 
atmosphere is forced to be hydrostatic by the stellar wind. 



side mass loss has been quenched, we have implicitly as- 
sumed that the planet resides outside of the region where 
the dipolar component of the star's magnetic field dom- 
inates. Inside this region — i.e., inside the star's Alfven 
radius — most magnetic field lines are not blown open by 
the stellar wind. Closed, poloidal field lines can encage 
both dayside and nightside winds, confining atmospheric 
escape to magnetic flux tubes in the vicinity of the 
planet's poles. The detailed structure of the solar mag- 
netosphere h as not been conclusive l y established. Some 
models (e.g.. iMcKenzie et al.l [19971 iMarsch et al.l |2003t) 
place the solar Alfven radius near 10 i?©, near the orbit 
of a hot Jupiter. However, th e three-dimensional model 
of iBanaszkiewicz et al] ()1998l their "DQCS" model with 
Q = 1.5) states that at lOi?0, field lines are primarily ra- 
dial, and by implication, the Alfven surface lies at smaller 
radius. 

What about for hot Jupiters orbiting T Tauri stars? 
To what extent are their dayside outflows squashed by 
the highly pressurized winds emitted by young active 
stars? If magnetospheric truncation of T Tauri accre- 
tion disk s sets the final o rbits of inwardly migrating hot 
Jupiters (jLin et al.lll996f ) , these planets likely also reside 
near the Alfven radii of their host star magnetospheres, 
where the stellar ram and magnetic pressures are com- 
parable. This conclusion is compatible with observations 
of mass loss rat es from T Tauri stars in the range 10 ~® 
to lO-'^Mg/yr (|Edwards et al.lll987[) and T Tauri sur- 
face fields that are 10^ tim es stronger than their m ain- 
sequence counterparts (e.g.. lJohns-Krull et ani999f ). We 
therefore estimate that P*_tot ~ 10 x (10"^)^ pbar ~ 10 
/zbar. Let us compare this pressure to the pressures char- 
acterizing planetary outflows. Figure fT3l shows a family 
of breeze solutions for a T Tauri-like UV flux of Fuv = 
5 X 10^ ergcm"^ s""'^. It shows that max(Ptot(?'s)) ~ 4 
nanobar <C P*,tot- Thus, it seems safe to conclude that T 
Tauri stellar winds completely stifie dayside winds from 
hot Jupiters. 

It is not clear whether hot Jupiters around T Tauri 
stars reside inside or outside their host stars' Alfven 
radii. If stellar rotation rates are locked to disk 
rotation rates at Alf ven radii ("disk locking"; e.g., 
iHerbst fc Mundtll200 Ell. this question reduces to whether 
the planet orbits inside or outside the corotation cir- 
cle. The current m easured rotation period of HD 
209458 of ~12.4 days (jWinn et al.ll2005l and references 
therein) yields a corotation radius of '^0.1 AU. How- 
ever, T Tauri stars typically ro tate faster than their 
main-sequence counterparts ( e.g., iHartmann et al.lll986l : 
iJohns-Krull fc GaffordI 120021 ). Furthermore, observa- 
tions of T Tauri magnetic fields do not consistently con- 
firm models in wh ich the closed field extends to the 
corotation radius (jJardine et al.l l2008l and references 
therein). 

Given these uncertainties, we acknowledge two possi- 
bilities. Winds from giant planets located well inside the 
magnetospheres of their T Tauri host stars will be con- 
fined on both day and night sides by magnetic pressure, 
blowing only along polar flux tubes. Hot Jupiters located 
in open fleld line regions experience quenching only on 
the day side — they lose mass strictly by night. 

3.5. Radiation Pressure 
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Fig. 12. — Mass loss rates (assuming 47r symmetry) for our breeze solutions. Mass loss rates become arbitrarily small as the breeze 
parameter [3 — which sets the maximum breeze velocity at Tb — decreases. I nter estingly, the total pressure Ptot('"s) asymptotes to a maximum 
value of ~40 pbar as ever slower breezes are considered (see also Figure [TTJ. The existence of a maximum pressure suggests that if the 
stellar wind exerts an external pressure in excess of ~40 pbar (for our standard parameters and boundary conditions), it will shut down 
the dayside planetary outflow completely. 

two orders of magnitude lower than the mass loss rates 
we have computed for photoionization-heated, hydrody- 
namic winds. We conclude that radiation pressure is not 
a significant driver of planetary outflows when compared 

a gainst photoionizatio n. 

IVidal-Madiar et al.l (|2003l ) claim a mass loss rate of 
j^Qio gg-i i3asg(j line-driven acceleration to large blue- 
shifted velocities of gas that has been raised, presumably 
by other means, to the altitude of the Roche lobe. Their 
claim is based in part on an assumed hydrogen density of 
no '--^ 2 X 10^ cm~^ at a distance of r = i?Rocho/2. But the 
corresponding column noi?R,ocho/2 ~ 4 x 10^^ ^ A^Lya 
would be optically thick to Lya photons (xLya ~ 20). 



Can stellar radiation pressure drive substantial plan- 
etary mass loss? Neutral H atoms absorbing stellar 
Lya photons feel a radiation pressure force i^LyafLya/c, 
where c is the speed of light. We use the solar Lya 
flux scaled to a = .05 AU, FLya 2.4 x lO^erg/cm^/s 
(jWoods et al.l l200"ol ). If the line were only thermally 
broadened with a velocity of ~10kms~"'^, the absorp- 
tion cross section at line center would be CLyQ ~ 6 x 
10^ ^''cm^. To account for extra Doppler broadening 
from a range of bulk velocities extending up to Vorh, 
we lower the thermally broadened cross section by an- 
other factor of 10: ULya ~ 6 x 10~^^cm^. Putting it 
all together, we find that the radiation pressure force 
is comparable to that of ste llar gravity, GM^mB _ /a? , i n 
agreement with the result of IVidal-Madiar et al.l (|2003D . 
But, radiation pressure of this magnitude would require 
of order an orbital period (days) to accelerate H atoms to 
velocities of Uorb ^ lOOkms"^. The observations, in con- 
trast, require this acceleration to occur over a few hours, 
the time to travel ~10i?p, the size scale probed by the 
transit. Furthermore, atoms are subject to Lya radia- 
tion pressure for only several hours before they are pho- 
toionized. These considerations imply that a larger Lya 
flux than we have assumed would be required to accel- 
erate neutral hydrogen to ^^100 km/s (jHolmstrom et all 
[20081) . 

Nonetheless, we assume optimistically that radiation 
pressure can in fact accelerate H atoms to '^Worb and 
estimate the resultant mass loss rate. Radiation pres- 
sure can only act on gas that is optically thin to Lya 
photons. The column density to optical depth unity is 
N-Lya = l/ffLyQ ~ 2 X lO^^cm"^. Hydrogen is acceler- 
ated off the planet's limb; the area of the annulus pre- 
sented by hydrogen within the planet's Roche lobe is 
~-^Rocho- This hydrogen is removed every i?Rochc/worb 
time. Putting it all together, we estimate a mass loss 
rate due to radiation pressure of 

n,T A^Lya^Rocho"^H Tn^„^-1 (01 \ 

jWrad. press. 7Z ] ~ iU gS , (Zi) 

^Roche/'yorb 



4. SUMMARY AND COMPARISON WITH OBSERVATIONS 

We have presented a simple model of atmospheric es- 
cape from hot Jupiters, as driven by photoionization 
heating by stellar UV radiation. To calculate the steady- 
state structure of the planetary wind, we employed a re- 
laxation code to solve the equations of ionization balance 
and of mass, momentum, and energy conservation. We 
imposed two-point boundary conditions that allowed us 
to find the unique wind solution that transitions from 
subsonic to supersonic velocities. Tidal gravity is impor- 
tant to include in our momentum equation because it 
alters the entire velocity structure of the wind, generat- 
ing higher velocities and helping to lower the wind's sonic 
point to within the Roche lobe, at least for the substel- 
lar streamline. Photoionization heating is balanced by a 
combination of PdV work and cooling by Lyman-a radi- 
ation. The latter serves as a thermostat, limiting plane- 
tary wind temperatures to ~10^ K over four decades in 
incident UV flux. Notably, conductive transport of en- 
ergy is not important, by contrast to the thermospheres 
of terrestrial planets. 

Our assumption that mass loss from hot Jupiters takes 
the form of hydrodynamic winds rather than local Jeans 
escape is self-consistent; in our solution, escaping gas 
remains collisional at the sonic point. Showing that 
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Fig. 13. — Same as Figure [TT1 but for the high flux case where 
i^uv = 5 X 10^ erg cm~^ s~^ . Over all wind and breeze solutions, 
the reference pressure Ptot(''s) does not exceed ~4000 pbar (middle 
panel). Since T Tauri stellar wind pressures are about three orders 
of magnitude greater than this (see text), we expect the daysides 
of hot Jupiters orbiting pre-main-sequence stars to emit neither 
winds nor breezes, but to have radially hydrostatic atmospheres. 
Their nightsides, however, may well be losing mass at a rate of 
order lO^^gs"^, since we expect horizontal winds to carry hot 
photoionized plasma from the daysido to the nightside. 

this condition is satisfied at the sonic point is sufficient 
because the flow inside the sonic point is denser and 
hence still more collisional, while the flow outside the 
sonic point is supersonic and so cannot influence the flow 
inside — in a supersonic flow, one cannot shout upstream 
and be heard. We agree with iTian et al.l (|2005l ) and 
iGarcfa Mufiod (|2007f) that the flow is never in a regime 
where Jeans escape considerations are relevant (see ij3.ip . 

We find that a planet similar to the transiting hot 
Jupiter HD 209458b loses mass at a maximum rate of 
^lO^^gs"^ for a UV flux characteristic of low to mod- 
erate solar activity. For a UV flux near the Lyman edge 
that is ^2-3 times larger, as obtains during solar max- 
imum (jWoods et all 120051: iLean et al.l 120031 and refer- 
ences therein), the mass loss rate is ~2 x 10^" gs""'^ (we 
have shown that at these fluxes, mass loss is energy- 
limited and scales nearly linearly with UV flux). This 
maximum mass loss rate is tw o orders of magnitude lower 
than the lO^^g/s quoted bvlLam mer et al.l (l2003l) . and 
inconsistent with the claim by Bar affe et al.l ()2004D that 
planets lighter than 1.5 Jupiter masses at 0.05 AU from 
their stars e vapor ate entirely in < 5 Gyr. We agree with 
iTian et all ()2005[ ) (who idealize the wind as neutral). 



lYelld ()2006[ ). and lGarcfa Muiiozl (|2007D that the current 
mass loss rate from HD 209458b causes the planet to lose 
at most 1% of its mass over its 5 Gyr age. 

What explains the factor of 100 discrepancy between 
our maximum mass loss rate s an d the mass lo ss rates 
given bv Lammer et all (l2003f) andlBaraffe et all (12004) ? 
These authors, inspired by IWatson et al.l ( T981h who 
model mass loss from the highly conductive atmo- 
spheres of terrestrial planets, posit that outflows from 
hot Jupiters are energy-limited. Energy-limited flows 
are those for which a fixed fraction of the stellar UV 
radiation incident upon a planet's surface goes towards 
driving gas out of the planet's gravitational well. As we 
have shown, conduction is not important in the winds 
emitted by ho t Jupiters, and therefore the details of the 
calculation bv I Watson et al.l ()1981h are not transferable. 
Nonetheless, we find that for F\jy < 10'* erg cm~^ s~^, 
as obtains for hot Jupiters orbiting main-sequence so- 
lar analogs, hot Jupiter winds are indeed nearly energy- 
limited. The energy-limited mass loss rate can be written 
as 

. efuv X nrl 
Mlim ~ „ I (22) 



(jWatson et al.l ll98lL see also our Equation [19]) , where 
gas is bound to the planet below radius tq and the bulk 
of incoming UV radiation is absorbed at ri. The dif- 
ference between our calculated M and that derived by 
iLammer et all (|2003[ ) arises from two factors. First, we 
include a heating efficiency e < 1, which must account at 
least for the energy lost to photoionizing atoms. In our 
simple model, e w 0.3 (recall Equation 2]). Second, and 
more significantly, our calculation places tq and ri closer 
to l.li?p (see beginning of i|2l), the location of the r = 1 
surface to photoionization. ILammer et al.l (|2003f) take 
instead tq k. ri k ?>Ru by ap plying detailed formulae 
derived bv I Watson et al.l (|1981h . These formulae are not 
appropriately applied to the photoionized upper atmo- 
spheres of hot Jupiters, because in thes e environments 
conductive cooling is not significant. iGarcia Mufioj 
()2007| ) reaches essentially the same conclusion; see his 
section 3.5. 

For high UV fluxes > 10^ ergcm"^ s~*, like those in- 
cident upon hot Jupiters orbiting active T Tauri stars, 
mass loss ceases to be energy-limited. Most of a fast 
photoelectron's energy is lost to collisionally excited Lya 
radiation. By contrast to the case at low Fyv where 
photoionizations are balanced by gas advection and the 
neutral gas fraction remains of order unity at altitude, 
at high i^uv photoionizations are balanced by radiative 
recombinations. The transition to a nearly completely 
ionized flow is sharp, and the overall wind structure is 
reminiscent of a classic expanding H II region. Whereas 
in the energy-limited regime M is expected to scale as 
i^UY, in the radiation/recombination- limited regime M is 

1 /2 

expected to scale as F^y , essentially because that is how 
the number density of ionized atoms scales in a plane- 
parallel Stromgren slab. Our detailed numerical model 
yields power-law indices of 0.9 and 0.6 at low and high 
UV flux, respectively. 

We have demonstrated that above the t — 1 sur- 
face to photoionization, our wind model is insensitive 
to our chosen boundary conditions. This helps to jus- 
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tify our neglect of hydrogen molecular chemistry. The 
uncertainty generated by this omission is embodied in 
our choice of the base radius rmin, which by definition 
is that radius where hydrogen is predominantly atomic. 
Without modelling the chemistry of H2, we cannot be 
sure where r,„in is located. But we have shown that 
we can estimate it to sufficient accuracy (about 10%; 
see the order-of-magnitude discussion at the beginning 
of [J2|) that our mass loss rates are uncertain by at most 
factors of a few (see Figure [53] of Appendix The 
bulk properties of our wind solutions ar e in good agree- 
ment with those of lGarcfa Munod (|2007l) a n d our ioniza- 
tion profiles ag ree well with thos e of lYelld (|2004l ). Both 
iGarcia Muii"o3 (|2007D and|Yelld ^2Q0^ model hydrogen 
and helium mo lecular chemistry in hot Jupiter winds, 
and lGarcfa Muh oz (2007il includ e s con tributions from D, 
C, N, O, and CH. I Garcia Muiiozl ()2007D finds that if met- 
als are present in the wind in solar abundances, they can 
increase the mass loss rate by factors of a few by reduc- 
ing the effectiveness of cooling at depth (below the 
T = 1 surface to photoionization). Again, these details 
can be co nsidered bur ied i n our paramete r rmin . 

Neither IMS (|2004D nor lGarcia Munod (|2007t l consid- 
ered the effects of Lyman-a line cooling — which we have 
demonstrated is important for large UV fluxes — or of line 
cooling from metals such as O II, O III, and N II. If clas- 
sic H II regions are any guide, metal line cooling could 
exceed Lya cooling, lowering hot Jupiter wind temper- 
atures by up to a factor of 2 (jOsterbrocldll989f ). Metal 
abundances in the upper thermosphere are extremely un- 
certain, depending on unknown turbulent mixing coef- 
ficients ("eddy diffusivities" ) . In any case, we do not 
expect metal line cooling to change our results qualita- 
tively. It can only lower mass loss rates and wind veloc- 
ities somewhat, strengthening our main conclusions. 

It is often said that transonic winds are charac- 
terized by zero pressure at infinite distance, while 
subsonic breezes ha ve finite pressure at infinity (e.g., 
iGarcfa Muhod [20071 ). But in practice, when consider- 
ing how the outflow from the planet's dayside interacts 
with the outflow from its host star, whether the dayside 
emits a wind or a breeze does not require us to examine 
conditions at infinity. Rather, we evaluate conditions at 
the sonic point. Again, as with our criterion for hydro- 
dynamic escape, the sonic point serves as discriminant 
because once the flow achieves supersonic velocities past 
the sonic point, it cannot influence the flow inside the 
sonic point. If the total external pressure (ram, thermal, 
and magnetic) exerted by the host star's wind is less 
than the total pressure exerted by the transonic plane- 
tary wind at its sonic point, then the transonic wind solu- 
tion obtains. Otherwise, either the dayside emits a more 
gentle breeze or — if the external pressure exceeds some 
critical value — the dayside atmosphere does not escape at 
all but is forced to be in vertical hydrostatic equilibrium. 
To the extent that the host star of HD 209458b emits a 
highly variable wind like that of the Sun, we find that 
the stellar wind pressure is comparable to the planetary 
wind pressure at the sonic point (both are measured in 
tens of picobars). During violent coronal mass ejections, 
the stellar wind pressure may overwhelm the planetary 
wind pressure. Thus we expect that dayside winds from 
hot Jupiters orbiting main-sequence solar type stars will 
alternately turn on and off (with a duty cycle that might 



well favor the off state). When the dayside wind is off, 
we expect the nightside wind to turn on and pick up 
the slack — the absorbed UV energy now being used to 
power horizontal flows that carry photoionized plasma 
to the planet's nightside, which is shielded from the stel- 
lar wind and therefore immune to pressure confinement. 

Winds from T Tauri stars have magnetic and ram pres- 
sures that are some six orders of magnitude greater than 
their main-sequence counterparts. Were a hot Jupiter 
orbiting a T Tauri star to also emit a wind, the pressure 
at its sonic point would be greater than under main- 
sequence conditions, but only by about three orders of 
magnitude according to our model. Thus, a T Tauri 
stellar wind will entirely quash dayside outflows from or- 
biting hot Jupiters. Outflows from the nightside (if the 
planet is outside the stellar Alfven radius) or along polar 
field lines (if the planet is inside the Alfven radius) would 
carry away at most ~2 x 10^^ gs~^ x 10^ yr or ~6 x 10""* 
of the planetary mass. These quantitative considerations 
rule out speculations that mass loss is significant during 
the host star's youth (e.g., Baraffe et al. 2004). 

Returning to the main-sequence case, to what de- 
gree will a hot Jupiter wind absorb stellar Lya pho- 
tons and produce an observable transit signature? Fig- 
ure [14] shows, for standard model parameters, how the 
extinction varies with wavelength in the Lya line, as- 
suming the planet is in mid-transit and that its orbit 
is viewed edge on. The extinction is integrated across 
the entire stellar disk (to perform this integral, we ex- 
tend our wind model out to r = 10i?p w 2i?r{ocho)- 
Local line profile functions are Voigt profiles, which in- 
clude thermal and natural broadening. Out to Doppler 
shifts of ±30kms~^, comparable to bulk wind veloci- 
ties, the line is essentially black (the Lya emission pro- 
duced by the wind itself is negligible since densities are 
low enough that the n = 2 population of neutral hy- 
drogen is well below the LTE value). At Doppler shifts 
of iblOOkms"^, the extinction plummets to 2-3%. This 
result is at odds with the observational claim that HD 
209458b reduces the stellar Lya fiux by some 9-15% at 
±100 kms"^ while in transit (IVidal -Ma diar et all 120031 : 
lBen-Jaffell[2007l: iVidal-Madiar et al.ll2008fK 

To drive this point home, we present in Figure [TSl the 
predictions of our model for the Lya line during transit, 
computed by combining t he out -of-transit line spectrum 
from Vidal-Madiar et ahj (|2003l . taken from their Figure 
2) with our extinction curve (Figurc [T^. Over the wave- 
length intervals used bv .Vidal-Madiar et al.l (|2003D . our 
model produces a flux decrement of 2.9%. Shown for 
comparison is the observed in-transit line spectrum from 
IVidal-Madiar et"al] (|2003f ). Agreement between model 
and observation is poor: the model spectrum is hardly 
absorbed at large velocities, by contrast to the observed 
spectrum. 

What could explain this discrepancy between our 
model and the observations? Three possibilities present 
themselves: (1) Our model underestimates the column of 
neutral hydrogen generated by the wind during transit 
by a factor of ~3-5; (2) physics that is missing from our 
model generates a population of neutral hydrogen mov- 
ing at velocities larger than the bulk velocity of the wind; 
or (3) the observed flux decrement at ±100 km/s is due 
to some combination of geocoronal emission and intrin- 
sic stellar variability (see also Ben-Jaffel 2007): processes 
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Fig. 14. — Maximum fraction of stellar Lya radiation obscured by 
absorption in tlic planetary wind, plotted as a function of Doppler 
shifted velocity from line center. We apply our wind solution for the 
substellar streamline over all 4tt steradians of the planet's surface. 
We assume that the planet is located at mid-eclipse, and that its 
orbital inclination relative to the stellar equator is zero. For each 
projected distance from the planet's center ranging from Rp to 
8.2_Rp = we calculate a line-of-sight optical depth, integrating 
up to a maximum altitude of lORp. We include both natural and 
thermal line broadening. The ordinate is the fraction of the stel- 
lar disk that is obscured, accounting only for e~'^ absorption and 
neglecting scattering and wind self-emission. All of our assump- 
tions maximize the obscuration. Absorption drops sharply beyond 
a few tens of km/s to a constant of 1.5% (horizontal dashed line), 
the decrement in the visible continuum. The velocity scale of ~20 
km/s is set by thermal broadening (~10 km/s at 10'* K), together 
with the bulk velocity (10—30 km/s at Mach numbers of a few). If 
the wind is emitted only from the dayside (nightside), our extinc- 
tion curve will be valid only at wavelengths redward (blueward) of 
line center. Absorption by the planetary wind does not account 
for the claimed ~9-15% decrement in flux from HD 209458b at 
Doppler equivalent velocities near ±100 km/s, as integrated over 
the intervals bet ween the two dotted lines and between the two 
dot-dashed Unes l|Vidal-Madiar et al.|[2003 : iBen-Jafll'el .20071 ). We 
further emphasize this inconsistency in Figure [TSl 

that, if observed over long enough times, should not be 
correlated with the orbital phase of the planet. We now 
comment on possibilities (1) and (2). 

Our wind model generates a Lya flux decrement of 
~3% at Doppler-equivalent velocities of ±100km/s due 
to naturally broadened, Lorentzian line wings. If the col- 
umn of neutral hydrogen traversed by stellar Lya pho- 
tons along lines of sight located 5-10i?p from the planet 
were larger by a factor of 3-5, then absorption generated 
by the line wings might generate the 9-15% absorption 
observed at ±100km/s. Ben-Jaffel (2008) argues that 
this is the case using the wind profiles of Garcia-Munoz 
(2007), who calculates that the wind is '^30% neutral 
at ~5i?p. At these distances, we find a neutral frac- 
tion of ~13% for our base case and ~6% for parameters 
matching those used in Garcia-Munoz (2007), though in 
other respects our wind solutions largely agree. Our ion- 
ization profiles are in better agreement with the solu- 
tions of Yellc (2004), which do not generate sufficient 
absorption in the Lorentzian line wings to match obser- 
vations (Ben-Jaffel 2008). The reason for these differ- 
ences merits further attention. As previously noted, ID 
wind solutions such as those in this paper, Yelle (2004), 
and Garcia Munoz (2007) break down at planetoccntric 
distances greater than ^5-10i?p. Multidimensional cal- 
culations could yield a larger neutral column (J. Stone, 



Fig. 15. — Observed out-of-transit (dashed line) and in- 
transit (dotted line) spe ctra, reproduced from Figure 2 of 
IVidal-Madiar et al.l l l200a) . In the line "core" from -42 to -f32 
km/s, where interstellar absorption is strong, we set the flux 
to zero. Our theoretical in-transit spectrum (solid line) is com- 
puted by attenuating the observed o ut-o f-transit spectrum accord- 
ing to our extinction curve in Figure [141 Theory predicts substan- 
tially less absorption than is claimed to be observed. However, 
the s pectrally unresolved measurements of IVidal-Madiar et al.l 
II2004I ) . taken at a different epoch and in principle pertaining 
to much smaller Doppler equ i valent velocities than were impli- 
cated by IVidal-Madiar et al.l 120031) . may well be due to the 
wind. Integrated over the line profile outside the core, our model 
predicts a flux decremen t of 2.4%, close to that observed by 
IVidal-Madiar et al.l ll200l) . 

personal communication) . 

Alternatively, is there some qualitative physics that our 
model is missing that could generate a population of neu- 
tral H atoms moving at velocities of ^100 km/s? An ap- 
peal might be made to interaction of the planetary wind 
with the stellar wind: in the bow shock, neutral hydrogen 
from the planet mixes with fast moving stellar plasma 
and might be accelerate d to large blueshifte d velocities 
(see the simulations by IStevens et al.l |1992| of colliding 
stellar winds). But no similar argument can be made for 
the observ e d red shiftcd absorption — which, according to 
iBen-Jaffeil (|2007| ). sometimes appears stronger than the 
blueshifted absorption (see tr ansit B2 in his Figure 3b) . 
Appeals to radiation pressure (IVidal-Madiar et al.ll2003f l 
founder for the same reason. iGarcia Muiiozl ()2007| ) sug- 
gests turbulence in the planetary wind itself as a way 
to broaden the line. But to generate the large veloci- 
ties observed, the energy in such turbulence would need 
to exceed the thermal and bulk kinetic energies in the 
mean flow by a factor of '^100. Such energy requirements 
seem insurmountable. Coriolis forces can turn stream- 
lines that are initially in the plane of the sky into our line 
of sight, producing both redshiftcd and blueshifted gas. 
But the time required for gas to reach ~10i?p, the size 
scale probed by the transit measurements, is too short to 
produce line-of-sight velocities of '^100 km/s. In short, 
because the planetary wind stalls at the bow shock as it 
blows towards the star, stellar gravity cannot accelerate 
the flow to large redshiftcd velocities. 

More promising is the possibility that charge-exchange 
between stellar wind protons and H atoms in the plane- 
tary wind could generate high velocity neutrals. Holm- 
strom et al. (2008) find that H atoms will be accelerated 
by charge-exchange to velocities of ~100 km/s on the 
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assumption that the steUar wind interacts directly with 
the planetary magnetosphere at ~4i?p and that neutrals 
from the planet have been lifted to that height. In ne- 
glecting the planetary magnetic field, we have modeled 
planetary winds whose pressures exceed those of their 
magnetic fields. In this case, charge-exchange in the 
shock between the stellar and plan etary winds might like- 
wise accelerate H atoms (see, e.g.. lRavmond et al]|1998l 
for a discussion of charge- exchange in the bow shock of an 
infalling comet) . Hot Jupiter magnetic field strengths are 
uncertain but magnetospheres may compete with plan- 
etary winds for the dominant source of pressure at high 
altitude ( §3.4|) . Whether the stellar wind forms a shock 
with the planetary wind or with the planet's magneto- 
sphere may vary from system to system. 

We conclude that although UV radiation from main- 
sequence stars can drive hot Jupiter winds with mass loss 
rates of -lO^^ g/s, the source of the observed absorption 
detected at Doppler-equivalent velocities of ±100 km/s 
in HD 209458b remains uncertain, with several possible 
candidates. What does our model predict for sp e ctrally 
unresolved measurements? IVidal-Madiar etaU ()2004l ) 
collect light over all Doppler equivalent velocities and in- 
dicate a wavelength-integrated flux decrement of 5 ± 2%. 
We take the o ut-of-tra nsit line spectrum from Figure 2 
of Vidal-Madi ar et all ([2003) and reduce it according to 
the obscured fraction computed in our Figure 1141 Inte- 
grating over the range 1213. 7-1217. 7A, and excluding the 
line core between 1215.5 and 1215. 8A (velocities between 
-42 km/s and -1-32 km/s) inside of which interstellar ab- 
so rption practically extingu ishes the line (sec Figure 1 
of I Vidal-Madiar et al.l 120031 ). we compute a wavelength- 



integrated Lya flux decrement of 2.4% (compare to the 
flux decrement in the visible continuum, 1.5%). As a 
check, we ap ply our procedure to t he ob served in-transit 
spectrum of IVidal-Madiar et al.l (|2003D . finding a flux 
decrement of 5.3%, in good agreein ent with the 5.7% 
quoted bv lVidal-Madjar et all i|2004f l. While our model 
flux decrement of 2.4% is sensitive to our assumed outer 
cut-off radius (lOJJp), and is uncertain because our model 
breaks down there (we neglect the full stellar gravity field 
and Coriolis forces), it is never theless close enough to the 
observed decrement of 5 ± 2% (| Vidal-Madiar et al.|[200l 
that the spectrally unresolved measurements may well 
be probing a planetary outflow. We look to the Hub- 
ble Space Telescpe to reproduce this signature of a hot 
Jupiter wind after STIS is repaired or the Cosmic Origins 
Spectrograph (COS) is installed. 
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APPENDIX 

A. SENSITIVITY TO BOUNDARY CONDITIONS 

We demonstrate that our solution is insensitive to our choices for BC3 through BC6 ('! )2.2.2|) . That is, we show that 
our standard model represents a "quasi-unique" solution that hardly changes over large and physically realistic regions 
of input parameter space. We also show that while our solution does depend sensitively on rmin (which does not enter 
as a formal boundary condition but represents instead a global scale factor), that parameter is known sufhciently well 
that it introduces no more than a factor of 2 uncertainty in our determination of the mass loss rate. 

Regarding BC3, Figures [TBI and [TTl show that as long as the base density p(r„iin) is large enough that T(rinin) ^ 1, 
the solution is insensitive to p(rmin)- Our standard value p(rmin) = 4 x 10~^'^gcm~^ gives T{r^in) = 50 and so satisfies 
this requirement. 

For BC4, we set /+(r,„i„) to an arbitrary number <C 1. When T(ri„in) ^ 1, /+(''min) ^ 1 and our solution is not 
sensitive to its exact value (Figures [TH] and [TO]) . In addition, we have verified that our solution does not change if we 
replace BC4 with the requirement that photoionizations balance radiative recombinations at rmin- 

For BC5, we have chosen r(r,„in) = 1000 K for our standard value. As demonstrated in Figures [20l and [2T| our 
solution is insensitive to this choice as long as T(rniin) ^ 10^ K (the temperature such that t hermal velocities are 
comparable to the local escape velocity). All models of hot Jupiter atmospheres at depth fe.g.. iBurrows et al.ll2003l ) 
have this property. 

Regarding BC6, Figures [22l and [23l demonstrate that our solution is independent of our choice for T(rs) as long as 
it is <^ 1. This requirement is satisfied by the quasi-unique solution that appears when varying the other boundary 
conditions BC3-BC5. Thus, we have shown that our quasi-unique solution self-consistently demands T(rs) <C 1. The 
condition that T(rs) ^ 1 is also reasonable because the optical depth of material outside the planet's Roche lobe 
should be small. 

Finally, because we do not solve for the structure of the planet's atmosphere below rmin, we are not sure whether 
our standard value for r^i^ — Rp = 10^° cm corresponds correctly to our adopted base conditions. In other words, 
we cannot say with certainty whether our adopted base temperature of 1000 K and base density of 4 x 10""'^^ g cm~'^ 
are indeed reached at our assumed radius of r,„iii = 10^° cm. The value of the base radius is important since it helps 
determine the effective planetary cross section for absorption of stellar radiation and the local planetary gravity, both 
of which affect the mass loss rate — see Equation ([T9l) . Figure [24l shows how our results depend on rmin- The mass loss 
rate M is indeed very sensitive to rmin, changing by as much as a factor of 2 when rmin changes by just 20%. This is 
consistent with the A/c-Um oc scaling found in Equation (I19|) . However, our uncertainty in r^-^ is only about 10%. 
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16.-- Mass loss rate M and r{r^ii-^) as a function of 
boundary condition BC3: p(?'min)- The other boundary con- 
ditions BC4-BC6 are kept fixed at their standard values, and 
^"1X1111 = Rp- To calculate M, we apply our ID solution over 
a full in steradians. As long as we choose p{r^[-^-^) sufficiently 
large that T(rmin) 2> 1, M is insensitive to our choice. 



Fig. 17. — Wind velocity as a function of radius for different 
choices of BC3: p(rmin). For /9(rmin) = 4 X 10"^^ to 5 X 10"^^ 
g/cm'\ the solutions are indistinguishable from the solid line. 
For lower densities p{r^^^) = 2 X 10"^'* (top dashed line) and 
1 X 10^^^ (bottom dashed line), T(rniin) is not 2> 1 and the 
profiles are sensitive to BC3. 
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Fig. 18. — Mass loss rate M as a function of BC4: f+(r^i^). 
The other boundary conditions BC3, BC5, and BC6 are kept 



fixed at their standard values, and r„ 



Rri- To calculate 



M, we apply our ID solution over 4tt steradians. 
/+(ymin) < 1, M is insensitive to /+(rniin)- 
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Fig. 19. — Wind velocity as a function of radius for different 
choices of BC4: f+{r^in). For /+(rmin) = 10"^ to 10"\ the 
solutions are indistinguishable (solid line). For /-|-(»'min) = 1 
(dashed line), the solution differs only slightly. For that solu- 
tion, /_|_ drops to a small value at r just above r^i^- 
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Fig. 20. — Mass loss rate M as a function of BC5: T(rniin). 
Boundary conditions BC4 and BC6 are kept fixed at their stan- 
dard values, and Tmin — Hp. For BC3, we take p('^inin) — 
4 X 10-13 g/cm3 for T{rmin) > 800 K, 8 X lO'^^ g/cm^ 
for 800 K > r(rniin) > 100 K, and 3 X 10" g/cm^ for 
T{r^in) = 100 K. These adjustmen ts i n p(rmin) arc made to 
maintain T(ri„in) S> 1 (see Figures 1161 and I17p . To calculate 
M, we apply our ID solution over Ait steradians. As long as 
T(rinin) < 10"* K, M is insensitive to T(rinin)- 
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Fig. 21. — Wind velocity as a function of radius for different 
choices of BC5: T(rmin). For T(rinin) = 100 to 1000 K, the 
solutions are practically indistinguishable from one another. 
For higher temperatures T{r^in) = 2000 K (lower dashed line) 
and 6000 K (upper dashed line), the wind profiles change. 
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Fig. 22. — Mass loss rate Af as a function of BC6: T(rs). The 
other boundary conditions BC3-BC5 are kept fixed at their 
standard values, and r^i^ = Rp. To calculate M, we apply 
our ID solution over 47r steradians. Provided T(rs) <C 1 — 
which it self-consistently is for our quasi-unique solution — M 
is insensitive to T(rs). 



Fig. 23. — Wind velocity as a function of radius for different 
choices of BC6: T(rs). For T(rs) = 10"^ to 10"^, the solutions 
are indistinguishable from the solid line. For larger T(rs) = 0.1 
(top dashed line), 1.0 (middle dashed line), and 2.0 (bottom 
dashed line), the solution changes. But r(rs) < 1 is demanded 
by the quasi-unique solution that appears when varying the 
other boundary conditions BC3-BC5. 

We should choose a value for rmin that lies between the l-bar radius (say) and the radius where the bulk of the stellar 
UV photons are absorbed. From the estimates made at the beginning of fJH 1 < Tynin/Rp ^ l-l- Therefore in practice 
the uncertainty in M due to our uncertainty in rmin amounts to no more than a factor of 2. 

B. OTHER COOLING AND IONIZATION MECHANISMS 
The extra cooling mechanisms we considered and found to be negligible include collisional ionization 

A,ol = -1.3 X 10-21„+norl/2e-157809K/Tg^g^^^3g-l 

(|Blacklll981l) , recombination radiation (appropriate for Case A and thus an overestimate) 
(|Blacklll98H) . free-free emission 



-2.85 X 10"2^T1/2(5.914 - 0.5 InT + 0.0118iT^/^)nl erg cm" 



3s-i 



Aff = -1.426 X IQ-^VT'^^^^+ergcm-^s"^ 
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Fig. 24. — Mass loss rate M as a function of model parameter r^^i^. Boundary conditions BC3-BC6 arc held fixed at their standard 
values. To calculate M, we apply our ID solution over a full An steradians. Though the dependence of M on r^i^ is strong, the latter is 
known to within 10% (see text). 



where gtf ~ 1.3 is the Gaunt factor ()Spitzeij|l978l ). and conduction 

Acond = ^^(^r.-j (B4) 

where the thermal conductivity n = 4.45 x 10'*(T/10^ K)*^- erg cm"^ () Watson et al.|[T98ll ). In equations 

(jBip - (jB3|) . densities are in cm~'^ and temperatures are in K. Note that conduction can either cool or heat gas locally 
(and indeed both signs are observed — see Figures [3] and [6]) . 

We also considered how colli sional ioniza tions change the ionization balance. The coUisional ionization rate is given 
by dividing (IBip by -13.6 eV ()Blacklll98ll ). We found this contribution to be neghgible. 

C. ESCAPE OF LYa COOLING RADIATION 

To act as an effective coolant, Lya photons must be able to escape the wind. But the wind is optically thick to Lya 
photons. Radiative cooling is thwarted if before the photons escape by resonant scattering, they excite H atoms that 
subsequently undergo collisional de-excitation, converting photon energy back into heat. We show here that this is not 
a significant effect. 

Line photons escape by frequency redistribution: scattering into line wings where the Lya optical depth TLya is 
much reduced. The number of scatterings A^scat required for a photon to escape is given approximately by the inverse 
of the probability Pscat that an excited atom emits the photon at a frequency such that TLya < 1. We estimate this 
probability as 

Pscat ~ 2 / 4>{y)dv 
J fi 

where </) is the Voigt line profile function, accounting for natural and thermal broadening at T = 10^ K, and vi is the 
frequency for which TLya = Ij blueward of line center. The frequency vi is such that 



where 



0(line center) 
TLya (line center) 



, ^ (TLya (line center) - 
TLya (hue center) ~ — t ^ 3 x 10 t 

and CTLya (line center) = 6 x 10^^^ cm^ is the absorption cross section at line center. For a photoionization optical 
depth T '--^ 1, we find that A^scat ^ Psclt ~ 1 x 10^. This estimate neglects differential bulk velocities in the wind, which 
tend to decrease iVgcat, and the random 3D directions with which photons are scattered, which tends to increase iVgcat- 
Both effects are expected to be of order unity. 

While the photon is being scattered, it spends a time ^cx ^ -Ascat^2i^ "locked" inside an excited H atom (in the form 
of electron excitation energy), where A21 = 6.3 x 10^ s~^ is the Einstein A coefficient. By comparison, if we assume that 
collisional de-excitations are dominated by fast thermal electrons, the time required for an excited H atom to experience 
a collisional de-excitation is tcoi ~ ('T-+o'cnWth,c)~^, where (Ten ~ 10~^^cm^ is the electron-neutral de-excitation cross 
section and Uth,c ^ 400kms~^ is the electron thermal velocity. In Figure [21 we show that n+ < lO^cm"'^ for t < 1 for 
standard model parameters. Putting it all together, we find that ^cx/icoi ^ 6 x 10^'^. Thus the photon readily diffuses 
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out of the wind before becoming thermalized. Moreover, the photon spends such a smaU fraction of its time locked 
inside an excited H atom that outward advection of gas does not change this conclusion. 
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